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Abstract 



We solve the dynamics of an ensemble of interacting rotors coupled to two leads at 
different chemical potential letting a current flow through the system and driving it out 
of equilibriuixL. We show that at low temperature the coarsening phase persists under the 
'^ ■ voltage drop up to a critical value of the applied potential that depends on the charac- 

teristics of the electron reservoirs. We discuss the properties of the critical surface in the 
Q ! temperature, voltage, strength of quantum fluctuations, and coupling to the bath phase 

diagram. We analyze the coarsening regime finding, in particular, which features are 
CS| ■ essentially quantum mechanical and which are basically classical in nature. We demon- 

^ i strate that the system evolves via the growth of a coherence length with the same time 

dependence as in the classical limit, R{t) ~ i^" - the scalar curvature driven universality 
TJ" ■ class. We obtain the scaling function of the correlation function at late epochs in the 

coarsening regime and we prove that it coincides with the classical one once a prefactor 



in I that encodes the dependence on all the parameters is factorized. We derive a generic 

formula for the current flowing through the system and we show that, for this model, it 



rapidly approaches a constant that we compute. 



Contents 

1 Introduction 



The model 

2.1 System of disordered quantum rotors 

2.2 Reservoirs of electrons 

2.2.1 DOS witli a finite bandwidth 

2.2.2 DOS at low energy 

2.3 Coupling between the system and the reservoirs 



3 The dynamics 

3.1 Quench setup J9 

3.2 Schwinger-Keldysh formalism 

3.3 Macroscopic observables 

3.3.1 Keldysh Green's functions 

3.3.2 Self correlation 

3.3.3 Self linear response 

3.3.4 Keldysh rotation 

3.3.5 Bosonic FDT 



5.1 Average over disorder 

5.2 Schwinger-Dyson equations 

5.3 Quantum non-equilibrium steady state (QNESS) phase 

5.4 Critical manifold 

5.4.1 Equation for criticality 

5.4.2 Critical points on the F = plane 

5.4.3 Quantum critical point 

5.4.4 Summary of the phase diagram 

5.5 Coarsening phase 

5.5.1 Numerical solution 

5.5.2 Mapping to Langevin dynamics 

5.5.3 Long-time dynamics 

5.5.4 Linear response 

5.5.5 The Lagrange multiplier 



The influence of the fermion baths ll4 

4.1 Self-energy 

4.2 Some limits 

4.2.1 Zero drive 

4.2.2 Low frequency 

4.2.3 ep much larger than all other energy scales 

4.2.4 Zero temperature 

4.2.5 Some specific reservoirs 



14 
15 
15 
16 

m 

17 
U 



Results h8 



19 
20 
21 
21 
21 
23 

25 
26 
26 
26 
27 
29 
33 
33 



6 The current 35 



B.l Keldysh rotation 

B.2 Symmetry properties under t -^ t' 

B.3 Free fermions 

B.3.1 Single free fermion .... 
B.3. 2 Collection of free fermions 
B.3. 3 Fourier transforms .... 



D.l Derivation within the Schwinger- Keldysh formalism 
D.2 FDT check 



E.l Quadratic effective action 
E.2 Saddle-point evaluation . . 
E.3 Schwinger-Dyson equations 



7 Conclusions and discussion 38 



8 Acknowledgments |38 



A Conventions |39 

A.l Fourier transform 39 

A. 2 Heisenberg representation 39 

A. 3 Time-ordering operator 40 

A. 4 Green's functions 40 



B Fermionic bath kl 



41 
41 
42 
42 
42 
43 



Fluctuation-Dissipation Theorem |43 



D Computing the self-energy |44 



44 
47 



E Dynamics |48 



48 
49 
50 



1 Introduction 

Quantum mechanics determines the behavior of physical systems at atomic and subatomic 
scales. The search for quantum effects at macroscopic scales started soon after the development 
of quantum mechanics. A number of quantum manifestations at such scales have been found 
including quantum tunneling of the phase in Josephson junctions pLj or resonant tunneling of 
magnetization in spin cluster systems [2]. 

Dynamic issues in isolated quantum many-body systems are the focus of active research. 
Some of the problems that are currently being studied theoretically are: the time evolution of 
the entropy of entanglement in spin systems [3] , the nature of non-equilibrium steady states in 
small quantum systems driven out of equilibrium |H |5] due to their relevance for nano-devices, 
quantum annealing techniques [6], and the density of defects left over after a gradual change 
in a parameter |7j . The influence of an environment on the dynamics of quantum systems was 
also dealt with in a number of cases such as the spin-boson model [1] , disordered spin chains 
coupled to bosonic baths [8], or an electronic ring coupled to leads and further driven by a 
time-dependent field [9l [10] . 

Once the interest is set upon macroscopic systems, the question as to whether these undergo 
phase transitions naturally arises. The theory of equilibrium classical and quantum phase 
transitions is well developed. Non- equilibrium phase transitions in which quantum fluctuations 
can be neglected are also quite well understood. These are realized when a system is forced in a 
non equilibrium steady state (by a shear rate, an external current flowing through it, etc.) pT| 
[T2l|T3l|ll] or when it just fails to relax {e.g. after a quench) and displays aging phenomena [151 
[16] . In contrast, the effect of a drive on a macroscopic system close to a quantum phase 
transition is a rather unexplored subject. Some works have focused on non-linear transport 
properties close to an (equilibrium) quantum phase transition (TTJ [HI [19| . Others have studied 
how the critical properties are affected by non-equilibrium drives [20l [2T| [22] . However, a 
global understanding of phase transitions in the control parameter space T, V, F, with T 
the temperature, V the driving strength, and T the strength of quantum fluctuations, is still 
lacking. Furthermore, to the best of our knowledge, the issue of the relaxation toward the 
quantum non-equilibrium steady state (QNESS) has not been addressed in the past. 

In this paper we extend our study of driven quantum phase transitions and coarsening 
phenomena started in [23] . We study a class of analytically tractable models, systems of 
M-component N quantum rotors that encompass an infinite range spin-glass and its three 
dimensional pure counterpart modeling coarsening phenomena. As discussed in [2l| models of 
quantum rotors are non-trivial but still relatively simple and provide coarse-grained descriptions 
of physical systems such as Bose-Hubbard models and double layer antiferromagnets. The 
system is coupled to two different external electron reservoirs that lead to a current flowing 
through it and driving it out of equilibrium. (For a two-dimensional model the current flows 
perpendicular to it, see the sketch in Fig. 1 of [20].) In the simplest setting [20] each rotor 
is coupled to independent reservoirs; more realistic couplings are discussed in [22j. Using the 
Schwinger-Keldysh formalism [251 [2S] we obtain the complete out of equilibrium dynamics of 
these models in the large M limit. We show that at sufficiently low T, V, F, see Fig. 1, the 
system never reaches a QNESS and ages with remarkable universal properties. We study the 
critical properties of the phase transitions, in particular in the vicinity of the (drive-induced) 
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Figure 1: Non-equilibrium phase diagram of the fully connected driven quantum rotor model 
with an infinite number of components. 

quantum out of equilibrium critical point l^c at F = 0, T = and the "usual" quantum critical 
point Fc at \^ = 0, T = 0. We analyze in detail the relaxation in the coarsening regime and 
uncover the scaling properties of correlation functions and linear response. We derive a general 
formula for the current flowing through the system under such a voltage drop and we analyze 
its dependence on the dynamics of the system. Some of these results were announced recently 
in [231. 



2 The model 



2.1 System of disordered quantum rotors 

The model we focus on is a quantum disordered system made of N M-component rotors in- 
teracting via random infinite-range couplings [2^. We consider a fully-connected (mean- field) 
model where there is no underlying geometry: each rotor is equivalently coupled to all the 
others. The Hamiltonian is given by 
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n, • n. 



(2.1) 



nf (/i = 1 . . . M) are the M components of the i-th rotor. The coordinates n^ constitute a 
complete set of commuting observables. The scalar product rij ■ n^ is given by Ylii=i ^t^^- The 
length of rotors is fixed to unity: rij ■ n.j = 1, V i = 1 . . . A^. The strengths Jjj's are taken from a 
Gaussian distribution with zero mean and variance J^ . J controls the strength of disorder. Lj 
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Figure 2: Density of states (DOS) of type A reservoirs, hq and /Zq + eV are the left and right Fermi levels, 
respectively. The left reservoir is half-filled. 
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r acts like a moment of inertia and controls the strength of quantum fluctuations; when 
fi?T / J — 7- the model approaches the classical Heisenberg fully-connected spin-glass. In the 
large M limit it is equivalent to the quantum fully-connected p = 2 (or Sherrington-Kirkpatrick) 
spherical spin-glass [2HI [29]. The classical mapping to ferromagnetic coarsening in the 0{Af) 
model with M —> oo [16] holds, as we shall show in Sect. I5.5.3[ for the quantum model as well. 



2.2 Reservoirs of electrons 

The system is coupled to two, 'left' (L) and 'right' (R), reservoirs of electrons. These inde- 
pendent reservoirs are both in equilibrium at inverse temperature Pl and Pr. The situation 
Pl 7^ Pr would create a heat flow from one reservoir to the other. We are interested in the 
simpler case in which Pi = Pf/^ = P = T^^ {ks = 1). An electric current is forced by imposing 
different chemical potentials, iii = A*o and fJ^R = fio + eV (where — e is the electric charge of 
one electron). eV is the strength of the drive. As eV/J — )■ 0, the effect of the reservoirs on the 
system approaches the one of an equilibrium bath at temperature T. The details of the reservoir 
Hamiltonians H^ and Hr are not important since only the electronic Green's functions matter 
in the small rotor-bath coupling we concentrate on. We consider the simple case in which left 
and right fermionic reservoirs have the same density of states (DOS) pi = Pr = P- Moreover, 
we focus on simple cases in which the shape of the DOS is controlled by only one typical energy 
scale ep. In the rest of this paper, we often consider the limit in which ep is much larger than 
all the other energy scales involved. In this limit the results become independent of the detailed 
functional form of the DOS. We also give some results for finite ep using the specific DOS that 
we introduce below. 
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Figure 3: An example of type A reservoir: the semi-circle density of states (half- filled) . 

2.2.1 DOS with a finite bandwidth 

We first consider regular DOS which have a finite typical width (finite bandwidth) controlled 
by ep and /xq is set around the maximum of the distribution. In the limit where ep is very 
large, they can be seen as almost fiat distributions. We call ecut the finite energy cut-off beyond 
which the DOS vanishes, p(|e| > e^ut) = 0. Since the DOS we consider have a single energy 
scale ep, ecut should scale with ep. Notice that a finite ecut constrains the voltage not to exceed 
eVmax = ^cut — fJ'O siucc the right reservoir is then completely filled and therefore it cannot 
accept more fermions. 

We call reservoir of type A a half-filleclj reservoir the DOS of which has a finite bandwidth 
controlled by ep and is symmetric and derivable in the vicinity of its maximum (see Fig. [2]). 
The simplest example of a type A reservoir is given by the semi-circular DOS (see Fig. [3]), 
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(2.3) 



that is sjTiimetric and centered around ep. Here ecut = 2ei?. We choose /io = cf so that the 
reservoirs are half-filled at zero drive {eV = 0). In this case, at T = 0, the voltage applied 
between both reservoirs cannot exceed eVmax = ^cut — yUo = ^f- 

Type B reservoirs have finite bandwidth but no energy cut-off: ecut = eVmax — ;■ C)0. A 
realization of these reservoirs is given by the following DOS [see Fig. 111(a)] 
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(2.4) 



where a ~ 0.97 is a numerical constant fixed by normalization. The maximum of this distri- 
bution is located at ep/\/2. This reservoir is half-filled for /ig ~ 0.95 ep. This distribution 
resembles the semi-circular one in the sense that they both start with a square root behavior, 
have a maximum, and a bandwidth of order ep. In contrast, the DOS in eq. (12. 4p is different 
from zero at all finite e and one can exploit this feature to apply strong voltages. 



-'^ Half- filled means that half the total number of available states are occupied: /^° de/3(e) 



i at T 
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Figure 4: Two examples of type B reservoirs, (a) The distribution pB vanishes asymptoticaUy. (b) The square 
root distribution diverges asymptotically. 

2.2.2 DOS at low energy 

In the previous examples (p^ and pb), we focused on values of po corresponding to high energy- 
states where the DOS is regular. We are also interested in studying cases where po is centered 
around low energy states. To analyze these cases, we focus on a DOS which reads [see Fig. Hl^b)]: 
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(2.5) 



This square root behavior is actually the one of the 3d free fermions reservoir. In this case 
ep is of the order of the hopping term for the free fermions. Since we shall only focus on the 
low energy states of the reservoir, we can neglect the non trivial high energy structure of the 
reservoir and take the DOS equal to zero for e > 2eF- 

For the 2d free fermions, the density of states is given by 



Pc2d[e) 



1 
2^ 



whereas for the Id free fermions, the density of states is given by 
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(2.6) 



(2.7) 



and, as for pcsd, we take these two densities of states to be equal to zero for e > 2eF- 

2.3 Coupling between the system and the reservoirs 

An electron hop from the L{R) reservoir to the R{L) reservoir is linearly coupled to each rotor 
component: 
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X] Xl 5Z XI ^^^' ^i ^'^Liki ^U' i^Rik'i' + L^ R] 



(2.8) 



i=l fj,=l k,k'=l l,l'=l 



where "ipl^^j^i is the l-th. component of an A^-component spinor operator that creates an additional 
fermion with energy hwk in the L reservoir associated to the i-th rotor, k labels the electron 
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energy inside the reservoirs, Ng is the total number of states in each reservoir, o"'^ are the 
generahzed Pauh matrices for SU(A^) of dimension A4 x A4 with A4^ — 1 = M. They are 
chosen to be normahzed such that Tia^a'^ = 5^^. Vkk' are the rotor-environment couphng 
parameters chosen to be constant: Vkw = ^c- Hsb is 0{MN) invariant. 

3 The dynamics 

3.1 Quench setup 

The system is initially prepared (at times t < 0) in such a way that its initial configuration 
(at time t = 0) is neither correlated with disorder (Jj/s) nor with the reservoirs. This can be 
realized, for instance, by coupling the system to an equilibrium bath at temperature Tq ^ J, F 
so that any correlation in the system is suppressed. At time t = the quench is performed 
by suddenly coupling the system to the L and R reservoirs. These are supposed to be "good 
reservoirs" in the sense that their properties are not affected by the state of the system. 

This setup generates non-equilibrium dynamics at times t > for multiple reasons. First 
of all, the rapid quenching procedure puts the system in a non-equilibrium initial condition 
with respect to its new environment. Moreover, the latter is not an equilibrium bath but a 
bias drive the role of which is to constantly destabilize the system. Finally, as a consequence 
of its disordered interactions, the system of rotors experiences intrinsic difficulties to reach 
equilibrium. Indeed, even if it were embedded within an equilibrium environment it would 
show a glassy phase [2S1 EI] in some parts of the phase diagram. 

Since system and reservoirs are decoupled at times t < 0, the initial density matrix of the 
whole system is given by 

N N 

g{t = 0) = gs{t = 0) (g) QLi (g) Qm ■ (3.1) 

i=i i=i 

QLi/Ri corresponds to the equilibrium density matrix of the L/R reservoir associated with the 
i-th rotor. The system of rotors being prepared at very high temperature, its initial density 
matrix is the identity in the rotors space: 

gs{t = 0)cxls. (3.2) 

All these density matrices are normalized to be of unit trace. The t > evolution of the whole 
system plus environment is encoded in 

g{t) = U{t,0) g{0) [t/(t,0)]t, (3.3) 

where the unitary evolution operator is given by U{t, 0) = Te""^ ■^o '^* ^(* ) with H = Hg + H^ + 
Hr + Hsb and T the time-ordering operator (see Appendix |A|). We analyze the non-equilibrium 
dynamics using the Schwinger-Keldysh formalism (see [22] for a modern review) that we briefly 
introduce in the following lines. 



3.2 Schwinger-Keldysh formalism 

The Suzuki-Trotter decomposition of the two unitary evolution operators that appear in 

Z = hm Tr f/(r, 0) ^(0) [U{t, 0)]"^ = 1 , (3.4) 

yields a path-integral involving two sets of fields with support on two different branches. The 
first ones are time-integrated on a forward branch from t = to +oo. In the following, these 
fields carry a + superscript. The other ones are time-integrated on a backward branch from 
-|-oo to and carry a — superscript. These two branches constitute the Keldysh contour C, see 
Fig. [5l The identity (13. 4p can now be expressed as a path integral, 

Z = y"l)[n±,V'^,V^^] e^^ (n+(0),^+(0)|f5(0)|n-(0),V-(0)) , (3.5) 

where we collected all the nf° fields into the notation n", and all the fermionic fields ip^^ and 



era 



their Grassmannian conjugates into i/'" and i/> (with a = ±). (n+(0), xp (0)|^(0)|n~(0), tjj~{0)) 
is the matrix element of the density matrix which has support at time t = only. The action 
5* is a functional of all these fields: 

/•oo 

5 = 5^ a/ dt£([r^^V'^V'"];t)• (3.6) 

a=± "^0 

The Lagrangian is given by £ = £5 + Csb + £l + ^^r with 

Cs{[nJ,t) = ^ E "^(^)' + ^ E '^^^ <W ■ ^K^) ' (3-7) 

csB{[n^,r,'^y,t) = v^^ E <"w fe'W ^i;' rmMt)+L^ r] . (3.8) 

ifikk'W 

C-L and Cr are the Lagrangians of the free fermions in the L and R reservoirs. The index 'c' at 
the bottom of the integral sign in eq. (13. 5 p is here to remind us that the integration is performed 
over fields satisfying the constraint that each rotor has a fixed unit length: n"(t) = 1 V a, z, t. 
The path-integral formalism gives a nice way to restore an unconstrained integration over all 
fields n" by the introduction of Lagrange multipliers z^: 

fv[n^]= fv[n'^]l[6{l-nm') (3-9) 

>^= -^ i,t 

= y'l?K,zTexp|^^"dta|^^^f(t)(l-<(t)^)j . (3.10) 

where we used the integral representation of the delta function (see Appendix |A]) and collected 
the new auxiliary real fields z^ into the notation z^. In terms of a Lagrangian, this gives rise 
to the new term 



10 



^ 00 



Figure 5: The Keldysh contour C goes from to +cx) and then back to 0. The Keldysh action involves forward 
fields (that live on the +branch of C) that are time-integrated from to +cx) and backward fields (that live on 
the —branch of C) and are time-integrated from -|-oo to 0. 

3.3 Macroscopic observables 

We are interested in the macroscopic dynamics of the rotors after an infinitely rapid quench and 
we wish to give an answer to the following questions (among others). Does the system reach 
a steady state? Does a steady state current establish? What are the long-time dynamics? 
We first obtain an effective generating functional for the rotors by expanding the system- 
drive interaction up to second order in the coupling, integrating away the fermionic degrees of 
freedom, and averaging over the disorder distribution. 

Introducing the external real fields rj'^At) that we collect in the notation 77" (t) (a = ±), the 
generating functional Z[ri^] reads 

Z[n^] = fv[n^,z^,iP^,iP^] ei^[-*'-*'V'±,v;*,r,±] (n+(0),V'^(0)|f?(0)|n-(0),V"(0)) , (3.12) 
where we introduced the source term 

S^S + hJ2Ut $^$^nr(t)r/r(t). (3.13) 

a=± i ^ 

The generating functional obeys the normalization property 2[t7^ = 0] = Z = 1 which is a 
fundamental feature of the Keldysh formahsm in this setup (see eq. (13. 4p and Sect. 15.11) . One 
has 



« W) - 



(3.14) 

T7±=0 



z 6vrit) 

where we introduced the notation 

(...)^ [v[n^,z^,i,^,^^] ■■■ e^^(n+(0),V'^(0)|f^(0)|n-(0),V-(0)). (3.15) 



Notice that one can distinguish this bracket notation from the quantum statistical average that 
we denote similarly by the occurrence of Keldysh indices inside the brackets. However, they 
coincide in the case of one time observables, e.g. 

{<it)) = {nTm , (3.16) 

with a = -|- or — equivalently if the observable is time-reversal invariant. 
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3.3.1 Keldysh Green's functions 

We introduce tlie two-time Green's functions G^jl^j^{t, t'), defined on the Keldysh contour (a, b 
±), as 



(-rw<(to) = -i ^'^f'^' 



ab 
7±=0 



ihG,f(t,t'). (3.17) 



Z 6vr{t)H{t' 

nf" being real fields, one has the following time-reversal property 

G.i(t,t') = Gi(t',t). (3.18) 

In the operator formalism, the Keldysh Green's functions read 

ihG.f^^it, t') = Tr [Tc n^{t, a) n]^{t' , b) g{0)] , (3.19) 

where n^^(t, a) denotes the Heisenberg representation of the operator nf at time t and on the 
a-branch of the Keldysh contour. Tc is the time-ordering operator acting with respect to the 
relative position of (t, a) and (t',6) on the Keldysh contour C (see Appendix El). 

We define the macroscopic Keldysh Green's functions by summing over the N rotors and 
each of their M components 

^ N M 

G^\t,t')^-Y,Y.Gui,{t.t')- (3.20) 

1=1 /i=l 



From the identity (I3.19p . one establishes two relations between the four Green's functions 

G++(i,i') = G-+{t,t')Q{t-t') + G+-{t,t')Q{t' -t) , 

G—{t,t') = G+-{t,t')Q{t~t') + G-+{t,t')Q{t' -t) , ^ ' 

leading to 



(3.22) 



Gr ~r Gr — Cj ~r Cj , 

G++it,t')-G—{t,t') = sign{t - t')[G-+it,t') - G+-{t,t')] . 

3.3.2 Self correlation 

We define the macroscopic two-time correlation as 

1 ^ 1 

= I [G+-(t,t') +G-+(t,t')] = Y [G-it,t') + G^^t,t')] . (3.24) 

It is symmetric in its time arguments C(t,t') = C{t',t). Given the constraint n(t) ■ n(t) = 1, it 
is one at equal times: G{t,t) = 1. The two-time correlation function is the simplest non-trivial 
quantity giving information on the dynamics of a system. In particular, a loss of its time 
translational invariance (TTI) is a signature of aging. 
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3.3.3 Self linear response 

The response at time t of the observable nf to an infinitesimal perturbation performed at a 
previous time t' on an observable /f linearly coupled to nf is defined as 



R^it,t')^'^<^'^^ 



(3.25) 

/r=o 



with the modified Hamiltonian 

H^H-ftn^. (3.26) 

Causality ensures that the response vanishes if t < t'. We define the macroscopic linear response 
as 

^ N M 

The functional derivative with respect to fi{t') in eq. (I3.25P can be written in terms of the 
source fields r]^ {t') since /f appears to play a similar role in the action functional: 

' -U.^-.^V (3.28) 



Therefore we obtain a Kubo relation, stating that the response can be expressed in terms of 
two-time Green's functions: 



N M 



«(m')^4ee^^^ '"'" 






T7±=0y 



= G"- (i, t') - G'"+(t, t') with a = + or - equivalently 

= i [G-(t,t') + G+-(t,t') - G++(t,t') - G-+(t,t')] 

= fG'+-(t,t')-G'-+(t,t')]0(i-i') , (3-29) 

where we made use of the relations (I3.2ip . 

Finally the four Keldysh Green's functions G"^^(t, t') can be re-expressed in terms of a couple 
of physical observables (namely correlation and response): 

mG''''{t, f) = C{t, t') - J [aR(t\ t) + hR{t, t')] . (3.30) 

3.3.4 Keldysh rotation 

The Keldysh rotation of the fields is a change of basis that simplifies the expressions of the 
physical observables such as the correlation C and the response R in terms of Green's functions. 
Moreover the connection with the Martin-Siggia-Rose generating functional in the classical limit 
is more straightforward in this representation [261 EI]- One introduces new fields as 

nn) = n/ — n • 



'■J ' 
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and the inversion relation 



< = n!na^nf . (3.32) 



We define the Green's functions of these new fields as ihC'^it^t') = 1/N X]i=i(^i('^) ' ^i(^')) 
with r,s = (1), (2). We have 

inG(") it, t') = Cit, t') , i/iG(i2) (t, t') = -iR{t, t') , 

inG(2i) (t, t') = -iR{t', t) , i;iG(22) (^^ ^z) = . ^'^■'^'^' 

The fact that G^^^) vanishes identically is very general and can be tracked back to be a conse- 
quence of causality. The unit length constraint imposed on the rotor coordinates, n"(t) ■n"(t) = 
1, becomes an orthogonality constraint between the fields in the new basis, nj (t) ■ n^^ (t) = 0, 
and a relation between their norms: n^^ (t) + ^Hj (t) =1. 

3.3.5 Bosonic FDT 

When the system of rotors is in equilibrium at a given temperature /3~^, the fluctuation- 
dissipation theorem (FDT) holds (in its bosonic version) giving an extra relation between the 
Green's functions. In Fourier space (see Appendix |S] for our Fourier conventions) it reads 

C{uj) = h coth {/3hio/2) Im R{u) . (3.34) 

For completeness, we derive this theorem in Appendix ID.2I 

4 The influence of the fermion baths 

4.1 Self-energy 

We treat the interactions with the environment in perturbation theory up to second order in 
the coupling. After the fermionic degrees of freedom are integrated out, the resulting effective 
action for the rotors acquires an extra term encoding the effects of the reservoirs. The detailed 
computation, given in Appendix ID.lt yields 

'S'cflF = Ss + Si^M + Sg^ , (4.1) 



with 



^4^][s«,«(^)] = ^M Yl /rdtdt'S-(t,t')X^n[(t) -11^(0, (4.2) 

rs=(l),(2) "^ "^° i=l 

and the four self-energy components 

Sf) = 2{n.:,r Re [Gf Gf - h^A (Gf Gf + G^,G^*)] ^ -Sf , (4.3) 

Sg^) = -2i(;^.)^ Re [Gf Gf + GM> iSi , (4.4) 

Sf) = 2i(na;J^ Re [G^Gf + GfGf] ^ -iS^ , (4.5) 

Sij'^ = 0. (4.6) 
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The fact that S^ vanishes identically is a consequence of causality. Similarly to what we have 
done in Sect. 13 .37^1 we renamed E^ , E^ and E^ into T,^, Eg and E;^. These real functions 
are usually referred to as the Keldysh, retarded and advanced components of the self-energy. 
G^, G^ and G-^ are the Keldysh, retarded and advanced Green's functions of the free electrons 
in the a-reservoir respectively (see Appendix IB. ip . Using their properties under time reversal 
(see Appendix IB. 21) . we establish 



^sir) 



^f( 



^R( 



-T 



(4.7) 



These relations reduce the number of independent self-energy components to two (namely E^ 
and E^). Plugging the expressions of the fermionic Green's functions given in Appendix lB.lt 
we obtain 



E^ 



^i( 



(r) 



li^cfii 



tanh(/3^^i^) tanh(/3^^^; 



^(^c)^(( 



tanh(/3 ) — tanh(/3 ) 



1 



sm 



cos 



Cl - e/?. 



h 



r )l)i 



(4.8) 






T)L)Reir). (4.9) 



The notation (( ■ ■ • )l)b. stands for fde^ de/j Pl{^l)pr{^r) ■ ■ ■ ■ The Fourier transforms read 



eKc) 



:7ih{hiJ,Y{{ 



(-L 



^^Uanh(/3^ 



PR. 



tanh(/3 ^ , ,,_ ^ 

X [5{nUJ - Elr) + 5{nUJ + eLR)])L)R 



(4.10) 



ReT.'^{uj) = -{nw,Y{{ 



tanh(/3 



^L- PL. 



tanh(/3 



Ci?, — Pr ^ 



pv 



^LR 



{hio) 



(-lr" 



;)l) 



R 



Im Ef (w) 



2^^(^c)^(( 



tanh(/3 



^L~ PL. 



tanh(/3 



e/?. — /^R - 



X [(5(^ - eLR) - 6{hw + eLR)])L)R , 



(4.11) 



where e^R = e^ — e^. Since E^(r) is a real and even function of r, E^(a;) is also a real and 
even function of cj. E^(r) being real, E|j(a;) is Hermitian: E^(a;) = E^(— a;)*. 

4.2 Some limits 



Expressions (I4.10p and (14.111) of the Keldysh and retarded self-energies are somehow cumber- 
some. We simplify them here in some physical limits. These expressions are heavily used in 
the rest of this work. 



4.2.1 Zero drive 

The L and R reservoirs constitute an equilibrium bath for the rotors as soon as they share the 
same temperature and the strength of the drive is set to zero {pl = Pr, eV = 0). In this case, 
the fluctuation-dissipation theorem applies to the bath, and gives an extra relation between the 
bath self-energy components. It reads 



^K 



,hw\ 



^R, 



B\^) = ^ coth ( I3-— j Im E^i^wj 



(4.12) 
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Ultimately the number of independent self-energy components reduces to one. We checked in 
Appendix ID.2I that the expressions (14.101) and (14. lip comply with the FDT in the equilibrium 
case. 

4.2.2 Low frequency 

Let us consider the low frequency limit (a; — t- 0), or long time-difference in real time, of the self- 
energy components of a generic non-equilibrium bath {eV ^ a priori). Parity considerations 
on T^^ and T,§ show that T,^{u) approaches S^(c<; = 0) which depends on T, eV and ep 
whereas Im Il^{u) oc u. The low frequency limit, which can also be seen as the classical limit 
{hu ^ T) of the quantum fluctuation-dissipation theorem in eq. (14.121) gives a way to express 
the temperature of an equilibrium bath as 

T=limi-^i^. (4.13) 

By analogy with the equilibrium case, we introduce for non-equilibrium situations 

T* = lim - ^^^"^^ (4 14) 

t.^0 2 d^\m Sg(u;) ■ ^ ' ' 

We expect that the effect of the reservoirs on the long time-difference dynamics of the rotors 
is the one of an equilibrium bath at temperature T* . 

4.2.3 tp much larger than all other energy scales 

The reservoirs act as an Ohmic bath in the limit in which ep is much larger than the tempera- 
ture, the drive and hw (eV, T,t\uj ^ ep). Equation (14. lip with Ae = e^^ — e^ reads 

ImSKw) = ^TrihiOcf fde' fdAepie')pie' -Ae)[5ihio-Ae)-5ihio + Ae)] 

tanh(/3^) - te^uHP ^^'' - ^'\- ^' - '^ )] . (4.15) 

In the limit hw <^ ep, we use p(e' ± hu) ~ p{e') and we derive 



ImJ:^^{co):^^n{huj,)'lde'p\e') 



tanh(/3^^ + ^-^°-^^) - t.nH/-^'^'^'^: 



The factor within the square brackets in the integrand is peaked at e' = po + eV. Hence we 
can approximate p^(e') ~ p'^{po) and then compute the remaining integral exactly to obtain 
an Ohmic (in the sense that it is proportional to u) behavior for the imaginary part of the 
retarded self-energy: 

Im S^(a;) ~ 2nh{fkOc)^p'^{po) w . (4.16) 

Interesting enough, this expression is independent of T and V. Similar calculations give 

^A-/ N ^/^ n9 9/ , e\^sinh(/3ey) — /iu;sinh(/3^) 

Sf a; ^ 2nh{nw,)'p'{po) tt^-^ ..n^\ ' ■ 4.17 

cosh(pe\/ ) — cosh.[prujj) 

In order to determine T*, we investigate the low frequency limit of S^(a;) given in eq. (14.170 . 
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Zero drive. For eV 't^T <^ ep, eq. (14.171) yields 

Sf (w) ~ 2nh'^{hLOcyp'^{fXo) ujcoth{l3nio/2) . 



(4.18) 



Equations (14.161) and (14.181) are linked through FDT. In the low frequency limit {hw, eV <^ 
T <t^ ep) it reads 

Sf (a;) ^ Anhihw^fp^ifio) T , (4.19) 

yielding T* = T as expected in this equilibrium situation. 



Finite drive. As soon as the drive is not negligible compared to temperature, in the low 
frequency regime {hw <^T <^ ep and eV -C ep) 



yielding 



E^(w) ~ 2'Kh{nuJc) p {pq) el/ coth (/3e1//2) 



T* = ^ coth {PeV/2) . 



An "FDT like" relation is verified in these limits 

Sf (w) = ;icoth(;ku/2T*)Im Sf(w) . 



(4.20) 
(4.21) 

(4.22) 



A similar interpretation of the effect of a two-leads bath in these limits on the dynamics of a 
single localized spin was given in [33] and |34j . 

Furthermore, in the low temperature limit {fiw ^ T ^ eV ^ ep) 



Sf (u;) ~ 2-Kh{f)w,f p\p^) \eV\ , 

yielding T* = \eV\/2. 

Finally in the zero temperature limit (0 = T ^ huj, eV <^ ep) 

Sf (c) = 2vrM7za;,)V^(/io) | j^| ,^ \J^ - \^y\ [ 



(4.23) 



(4.24) 



In the low frequency regime, we recover expression (I4.23p . In the zero temperature and zero 
drive limit (0 = T = eV ^ hw ^ ep) the Keldysh component of the bath self-energy reads 
T,^{u) = 27!-h{hUc)'^p'^{po) 1^1 that goes linearly to zero in the ^ — )■ limit. 



4.2.4 Zero temperature 

In the T = limit, we obtain for finite values of the other parameters (eV, hw, ep) 



Sf(c^) 



7rh{hwc 



Im J:^{u) = n{hujc)^ 



sign(eV^ -|- fko) / de p{e)p{e — hu) 

"'mo 

/•/j-o+sV—tiw 

+ sign(e\/ — hid) / de p{e)p{e + fku) 

liQ+eV+hu) pfio+eV—hu 

de p{e)p{e — fku) — de p{e)p{e + hw) 

MO "'mo 



(4.25) 
(4.26) 



17 



In the low frequency limit {0 = T <ti hu <t^ eV, ep) they yield 



Sf(w) ~ 2TTh{huc)^ sign{eV) de p^(e) , (4.27) 

ImE^iu)- 7rh{hw,f[p\fio)+p\fio + eV)]uj, (4.28) 

so that 

T*{T = 0) = sign(eV-) ^/^° 2/ V m • (4-29) 

4.2.5 Some specific reservoirs 

For the half-filled semi-circular DOS (type A), at zero drive and zero temperature, we establish 
the following analytical results at finite ep'- 

jnu:,\ \jl{Tep/h)-Sl[Tep/h) 

8 (nWc\^ Ji{TeF/h)Si{TtF/h) . , . 

B(r) , (4.31) 







Sf(- 






^""Bir. 


with Sg(r = 


= 0) = 


= 0, Sf (r = 



h\eF J {r/hy 

0) = ^{hjjJcY. Ji and 5*1 are the Bessel and the Struve functions 
of first kind and first order, respectively. From eqs. (I4.30p and (I4.3ip . we see that the temporal 
extent of both E^ and Eg is of order h/ep- In the limit in which ep is much larger than any 
other energy scale, a numerical analysis shows that this property holds for finite values of the 
temperature and the drive as well. As a way of summary, in Fig. [6] (a) we plot E^ as a function 
of rep for ep = 10 J, 100 J and at {T = J,V = 0) and (T = 0, ^ = J). In the case in which ep 
is finite, one can compute T* for the half-filled semi-circular DOS at zero temperature: 

V ) 2 1-1/2 {eV/epf \ \ max t y > 

In Fig. |6] (b) we give a numerical integration of Im Tj^iu) for the three types of reservoirs 
we introduced in Sect. 12.21 and in the case in which ep is the largest energy scale. This shows 
that the self-energy is indeed the one of an Ohmic bath. The fact that their Ohmic behavior 
is approximately valid until hw = ep supports the property that the temporal extent of the 
self-energies (in real time) is of the order of h/ep. 

5 Results 

In this section we present our results. We first complete the calculation of disorder averaged 
generating function and, from it, we derive Schwinger- Dyson equations for the two-time corre- 
lation and linear response valid for all values of the parameters. We next derive the dynamical 
phase diagram as a function of the temperature of the reservoirs (T), the strength of quantum 
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Figure 6: (Color online.) (a) S^ (for the half-filled semi-circle DOS) as a function of rep in the regime 
where tp is much larger than any other energy scale: for Ptp = 100 and (iep — 1000 at eV — and also for 
eV/2 — e_F/100 at T = 0. The three curves are indistinguishable. This shows that SH is indeed a function of 
rep in this regime and shows furthermore that eV/2 plays the same role as T. (b) Im I]]j(aj) is represented in 
a double logarithmic scale for the three following DOSs with Pep = ep/eV = 100: the half-filled semi-circle 
/OA(e), the half- filled type B with psi^) and the 3d free electrons DOS jOc3d(e)- The straight line above all is a 
guide to the eyes for a pure Ohmic (ex a;) behavior. The rapid decay above huj ^ ep \s a, signature of the energy 
cut-off, tcut c>c ef, of the DOS. 

fluctuations (F), the voltage {eV) and the coupling to the leads for which we introduce the new 
dimensionless parameter g = hwc/ep. We distinguish two phases separated by a second order 
phase transition. For high values of the temperature and/or strong drive and/or strong quan- 
tum fluctuations, we find a non-equilibrium steady state that approaches the usual paramagnet 
when eV — )■ 0. Whereas for low temperatures and/or low drive and/or quantum fluctuations 
we find a coarsening phase. 

5.1 Average over disorder 

At this stage, after tracing out all fermionic degrees of freedom, the effective action of our 
system is quadratic in the fields and reads 



i5e« = M 5: / dt a n.«(t) . r^it) + j=j: J, [4\t) . nf\t) + nf )(t) • n«(t) 
-\ j dt' Sf (t - t') nf (t) . nf (0 + i I dt' Si(t - f) nf (t) • r^\t'] 



a=± 



(5.1) 



Given that the initial condition for the rotors is taken to be uncorrelated with the disorder 
configuration (the Jjj's), neither the initial density matrix ^(0) nor the generating functional 
without sources (2[?7^ = 0] = 1) depend upon disorder. This property allows us to write 
dynamic equations by averaging over disorder the generating functional itself hence without 
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resorting to the use of replicas [3T]. As in other quantum systems with quenched disorder 
[29| |3T| [35| |37] , we are therefore interested in 



\id<i / 



(5.2) 



where P{Jij) is the Gaussian density distribution for the rotor couphngs with zero mean and 
variance J^. The disorder average over a random Gaussian potential can be readily done and 
the effective action of the system is quartic in the fields and reads 






h 



i=l 



(5.3) 



J'^M 



2N J dt' EK^W-f Wj K)(t')-nf (t') + nf)(t')-n«(t' 
i 

-\j<\t: Eg(( - (') nf'(() . nf ((') + xj<\t' Y.1(t - f) nf («) ■ «!'>(«' 



a=± 



5.2 Schwinger- Dyson equations 

In the large M limit, we show that the Lagrange multipliers are homogeneous, 

See Appendix [E] for a detailed computation. Moreover, introducing 



^K 



J^C + T. 



K 
B 1 



^R 



rR + T. 



R 

B 1 



(5.4) 



(5.5) 



we obtain the Schwinger- Dyson equations which fully determine the dynamics of the system: 

C{t,t') = [ dt"i:''{t,t")R{t',t")+ [ df^'^{t,f)C{t",t') , (5.6) 
Jo Jo 

R{t,t') = 6{t-t')+[dt"J:'^{t,t")R{t",t'), (5.7) 

Jt' 



'1 d^ 


+ 


zit) 


Tdt^ 


+ 


zit) 



z{t) = J dt" j:^{t,t")R{t,t") + j:^{t,t")C{t,t") - \^{t,t' ^ r) . 



(5.8) 



We remark that the expression for the response is decoupled from the self correlation apart 
from a residual coupling through the Lagrange multiplier. This is actually a consequence of 
two features of the model: the disordered potential is quadratic in the rotors and the coupling 
to the reservoirs is linear in the rotors. The "initial" conditions are given by 



C(t,t) = l, R{t,t) = Vt 



(5.9) 
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Moreover, integrating eqs. (15 .Op and (15. 7p over an infinitesimal interval around t' = t, one sees 
that the first derivative of the correlation is continuous at equal times 

lim dtC{t,t') = lim dtC{t,t') = , (5.10) 

whereas the one of the response function is discontinuous 

lim dtR{t,t') = r, lim dtR{t,t') = . (5.11) 

The structure of these equations is the same as the one in other out of equilibrium problems 
studied in [291 |3ll [351 |37l iO] . 

5.3 Quantum non-equilibrium steady state (QNESS) phase 

One expects that if the system is quenched into the high temperature phase, after a short 
transient it should relax toward a quantum non-equilibrium steady state (QNESS). The system 
of rotors cannot be in equilibrium since, for V^ 7^ 0, an electronic current is passing through 
it. Nevertheless the dynamics are still stationary (time translationally invariant). This implies 
that C{t, t') and R{t, t') are only functions of t — t'. Guided by a numerical analysis (see 
Sect. I5.5.5p . we make the assumption (that we later check to be consistent) that the quantity 
z{t) is a one-time observable that converges toward a finite value z°°. In this situation, one can 
Fourier transform the Schwinger-Dyson equations (15. 6 p and (15. 7p with respect to t — t' to find 

C{io) = J:''{io)\R{io)W (5.13) 

CM = JlgjI^^ Imi^M, (5.14) 

Using the fact that lim R^u) has to vanish, eq. (I5.12p imphes 

UJ—^OO 

R(cu) = ^ (-T-W + z'^- Sg(u;) + ^Ji-T-'u^ -|- z- - Sg(u;))' - A.P^ . (5.15) 

We note that in the cases in which the DOS of the reservoirs have an energy cut-off ecut, 

C{uj) = Im R{uj) = Ef (w) = Im S;|(u;) = ior hw > e^ut ■ (5.16) 

5.4 Critical manifold 

5.4.1 Equation for criticality 

Approaching the putative critical manifold from the disordered phase, see Fig. [H where after a 
short transient the system should be time translationally invariant, we look for a singularity in 
the Fourier transformed Schwinger-Dyson equations that would be the signature of the loss of 
time translational invariance and ultimately of a phase transition toward an out of equilibrium 
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behavior. Anticipating a second order phase transition scenario where the onset of criticahty 
is characterized by long- wavelength instabilities, we inspect these equations at a; = 0. 
The constraint that rotors have a unit length C{t,t) = 1 implies 

^ C{-) = i . (5.17) 

Q ZTT / 

and replacing C{u)) with its expression in eq. fl5.14p : 

"^^^l I,nflM = i. (5.18) 

Equation flS.lSp at w = reads 



Riu^ = 0) = ^ (^°° - ^Bi^ = 0) + V(^°°-Si(a; = 0)f-4J2^ 



(5.19) 



R{u] = 0) = /Q°°dr R{t) has to be real since R{t) is reajj. However, it is clear from eq. f l5.19p 
that z"^ = z^ = 2 J + S^((X' = 0) is a singular point (a minus sign would be incoherent with 
the approach in Sect. 15. 5p . This is the signature of the phase transition we were looking for. 
At criticality, 

Rico = 0)1^^^^^ = !/ J . (5.20) 

Concomitantly, the value of C{uj = 0) blows up. Inserting z'^ in eq. (I5.18p . we obtain the 
equation for the critical manifold, 

— ^%^Imi?(c.)|.o.= i. (5.21) 

27rImS^(w) V ;u. 2 ^ ' 

The parameters are the strength of quantum fluctuations F, the temperature T, the voltage 
applied between the two reservoirs V. We recall that J is the typical interaction between two 
rotors. The energy variation scale of the reservoirs is characterized by ep and huc quantifies the 
coupling strength of the rotors to their environment through the dimensionless small parameter 
g = hwc/ep. 

In the rest of this Section, we use eq. (I5.2ip to uncover the phase diagram of Fig. [TJ The 
critical surface is parametrized in the T, T V space by Tc, Tc, Vc [g is kept constant). We 
introduce the critical points % = Tc{T = V = 0),Vc = Vc{T = F = 0), f^ = Fc(T = V = 0). 
Anticipating the coming results, we introduce the dimensionless reduced parameters 6 = T/J, 
and V = eV/2J, 7 = (4/i/37r)^ T/J. In the plane ^ = 0, where the reservoirs act like an 
equilibrium bath, we recover the results in [21]. In the classical limit \^ = F = 0, we recover 
the ones in [30] . 

In the limit in which ep is much larger than any other energy scale, using eqs. (14.160 and 
(I4.17p . the equation for the critical surface reads 

do; 1 eVsinh(/3eV^) — ^sinh(/3^) of w ^ 

2^u cosh(/3eV^) - cosh(/3M ^^ ^(^)l^-= 2 • (^.22) 
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Figure 7: Study of the behavior of the 7 = critical Hne with the ratio ep/J for the half-fiUed semi-circle DOS. 
(a) The 7 = critical line dc{v) is given for four different values of the ratio ep/J. The analytical expression 
of the ep/J -^ 00 curve is given in eq. (|5.26p . For ep/J < 3/2 the critical point Vc is rejected to infinity, (b) 
Vc = Vc{d = 7 = 0) is plotted against ep/J. All these 7 = results are independent of the value g. 

5.4.2 Critical points on the F = plane 

Taking the F — )■ limit of expression f lS.lSp one has 



Im R{u') 



= 


it''- 


-(1- 


- a;'2)2 


for u' e [0, 72; 
for Lo' >V2, 



(5.23) 



where we introduced co' = u/V^JT. The expression of Im R{co) does not involve the reservoirs: 
the time scale of the rotors (controlled by F) totally decouples from the one of the reservoirs 
in such a way that the rotors only couple with the zero mode (the slowest) of the reservoirs. 
Using eq. (I5.2ip . we write the equation of the critical manifold in the F = plane 



lim 

r-5>o 




2F 



V2 



dcu' 



1-(1 



OJ 



i2\2 



Sf(V2jfw') 



27r ^ " ^" " ' Im T.^{^/2jfuj') 2 ' 
Using the definition (I4.14p of T*(T, eV) introduced in Sect. I4.2.2[ this simply reads 

T*{T,,eV,) = J. 



(5.24) 



(5.25) 



At eV = 0, for which the reservoirs constitute an equilibrium bath, the ratio S^/Im S|^ is 
given by the FDT and we find a temperature- induced classical critical point T^ = Tc{T = V = 
0) = J. In terms of the reduced temperature this reads 6c = 1- In the next two paragraphs we 
look at how this critical point is affected by a finite drive {eV 7^ 0). 



Infinite ep- We first consider the limit ei? — )■ 00, using the explicit expression (I4.2ip for T* 
one finds: 



TJeV) = — / arccoth ( — 
^ ' 2 / \eV 



(5.26) 



2s^B. 



I]^(a; — 0) is real for the same reason. 
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Figure 8: Phase diagram in terms of the reduced parameters analytically determined in the limit g 
(a) Critical line for V — 0. (b) Critical line for T = in the limit e^ — ?► oo. 



From this equation we find a drive-induced critical point at eVc/2 = J. In terms of tlie reduced 
voltage this reads Vc = 1- The departure from the classical critical temperature on the 7 = 
plane is quadratic: 6^ — 1 — (1/3) f^ for v <^ 1. Instead, on the zero-drive plane, v = 0, the 
critical line leaves 6c linearly: 6^ ^ 1 — (Stt^/IG) 7 for 7^1. More details on the critical line 
7c(t) at t; = are given in |27j. Close to Vc on the ^ = and 7 = planes the departure of the 
critical lines 7c('u) and Oc{v), respectively, are non-analytical and thus very steep [see Figs. [7] (a) 
and El (b)]. 



Finite ep. Let us now investigate the T = critical point Vc for finite values of ep. For 
our simple DOS depending on a unique parameter ep, Vc is controlled by ep/J. Plugging the 
expression fl4.29p for T*{T = 0) into the expression fl5.25p we obtain 



_ ^ 1 
sign(eV;, 



I 



Mo+eVc 

^J■o 



de'p2 



J p'^ifio) +p^(/io + eV;) 



1 . 



(5.27) 



The existence and the value of the solution Vc depend on the details of the DOS p(e). If the 
DOS has an energy cut-off ecut, the existence of a solution is guaranteed if the cut-off is larger 



than the solution e™" of 



dep\e) = Jp\fio) 



(5.2J 



MO 



For the type A half-filled semi-circle distribution (po = ^f, ^cut = Se^?), it turns out that 
eq. (15.271) admits a finite solution as soon as ep/ J > 3/2. For ep/J = 3/2, one finds eVc = 3/2 J 
{vc = 3/4) . For ep/J > 3/2, the finite solution Vc goes to one as one increases the ratio ep/J. 
For ep/J < 3/2 the critical point is rejected to infinity and the critical line in the F = plane 
converges to the asymptotic value 9c(v ::§>!) = 1/2 as ep/J — )■ 0. See Fig. [71 

For the distribution B, if po 7^ 0, the scenario is the same as for the semi-circle distribution 
there is a finite value of the ratio ep/J under which, the critical point Vc is rejected to infinity, 
and above which, Vc has a finite value that goes to 1 in the limit Ci? — > 00. If po = then Vc 
remains finite. 
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For the distribution of type C, eq. fl5.27p always admits a finite solution Vc independent of 
ep- For the distribution C3d, Vc = 1 regardless of /io, (-f and J . For the distribution C2d, we 
also get Vc = 1. For the distribution Cld, one can show that as long as /iq > 0, there is a finite 
Vc, function only of m = J/fiQi Vc = [exp {u + L(ue~")) — 1] /2m, where L(x) is the only solution 
of the equation Le^ = x that is analytic in 0. For ^^{ep — )■ oo) — )■ oo, we recover Vc = 1. 

5.4.3 Quantum critical point 

Weak coupling limit. We first consider the limit of the weak coupling to the reservoirs 
g ^ after the long-time limit such that the asymptotic regime has been established. It is 
actually in this (? — )• limit that the self-energy was computed (we expanded the total action 
up to second order in g) in Sect. HI (? = hwc/ep can be sent to zero by sending the coupling 
parameters to zero, but for our simple DOS, it can also be realized by sending ep to infinity. 
In equilibrium (l^ = 0) at T = 0, the FDT gives 



^\ , =hioiQ<nw < ecut ■ (5.29) 



Im S^fw 



By turning off the coupling to the reservoirs {g — )■ 0) in eq. (15.191) on has 



^ ^' ^ 1 for w'> 72, ^ ' 



where we introduced uj' = uj/\/2JV . Plugging eqs. (I5.29P and (I5.30p in the equation for the 
critical manifold ( ]5.2ip gives the quantum critical point 
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J if ecut > ~^J and no solution otherwise. (5.31) 



For type A reservoirs in the ep ^ oo limit, one can prove that the critical surface is parabolic 
close to the quantum critical point 7c, i.e., 7c ~ 1 — {16/3tt^)9'^ at 6^ ^ 1 and v = 0, and 
7c ~ 1 - (16/37r2) v^ for i; < 1 and ^ = 0. 

Finite coupling. When the coupling to the electronic reservoirs g is finite this quantum 
critical point (actually the whole critical surface) moves upward when increasing the coupling 
constant (see Fig. [9]). The coarsening phase is thus stabilized when increasing the coupling to 
the reservoirs. In the ej? — >■ oo limit, one has for g '^ 1 

7c^l + 2('^') {hwcfp'ifio) . (5.32) 

In the case of the type A half- filled semi-circle distribution this reads % ^ 1 + (9/2) g^. This 
is similar to what was found for other quantum spin models embedded in an Ohniic harmonic 
oscillator bath and is due to a spin- localization-like effect |29l |35]. This similitude is not 
surprising since we showed in Sect. 14. 2[ eq. (14.160 . that the mixed electronic reservoirs behave 
like an Ohmic bath in the ep ^ oo limit. 
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Figure 9: Numerical study of the evolution of the critical point 7c = 7c 
parameter g (here for ep/J = 10). 



= 0, w = 0) with the coupling 



rc(V=0) ~ Tc-T Tcir=0) ~ NA 

r4T=o) ~ NA Tc(y=o) ~ (Tc-r 



,1/2 



yc(T=o) 
K{r=o) 



-ryr- 



fT.-r)^/^ 



Table 1 : Behavior of the critical manifold close to the critical points for g 
critical point Vc = Vc{T = r = 0) the critical lines are non- analytical (NA). 



and ei? — >■ oo. Close to the 



5.4.4 Summary of the phase diagram 

Let us summarize the key features of the critical manifold in the case of a DOS with ep — > C)0. 
When the coupling to the reservoir g is set to zero, the values of three critical points {Tc, Tc, 
eVc) are only controlled by J that measures the disorder strength. Figure [1] gathers all the 
(yf — 7- results in the T,r,V space. The increase in either the thermal or quantum fluctuations, 
by raising F or the temperature T, respectively, leads to the destabilization of the coarsening 
phase. The same occurs for an increase in the bias voltage V. The summary of the behavior of 
the critical manifold close to the critical points Tc , Tc and Vc is given in Table [H Furthermore, 
an increase in the rotors-reservoirs coupling g pulls the quantum critical point Tc upward (as 
indicated in Fig. [1] by a vertical arrow) enlarging the low temperature phase. 

5.5 Coarsening phase 

We study the dynamics in the low T, weak F, weak V region of the phase diagram by solving 
the Schwinger-Keldysh equations in two ways: with an exact numerical approach and using 
analytic approximation in the long-time dynamics. We prove that in this region of the phase 
diagram there is coarsening and that the aging dynamics that occur are universal and equivalent 
to the ones of the classical (and undriven) limit of our model {a.k.a. the p = 2 spherical model 
with quenched disorder). 



5.5.1 Numerical solution 



Our numerical analysis consists in solving the Schwinger-Dyson equations (15. 6p . (15. 7p and (15. 8p 
after a quench into the low temperature, weak quantumness, weak drive phase. Thanks to their 
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causal structure, the equations on C, R and z can be integrated step by step in time, with a 
Runge-Kutta method. Apart from arbitrarily small numerical errors, this approach is exact. 

We concentrate on reservoirs at temperature T that have a type A semi-circle DOS (both L 
and R reservoirs). L reservoirs are kept half-filled while a voltage V is applied between L and 
R reservoirs, ep is chosen to be the largest energy scale. Typically, we consider the following 
values for the parameters: T ~ F ~ eV ~ O.IJ and ep ~ 10 J. 

The analysis shows (analytical arguments are given in Sect. I5.5.3P that the dynamics after 
the quench below the critical surface do not reach a QNESS. There is a separation of two- 
time scales typical of aging phenomena [l6] . The data in Figs. [T01IT2] were obtained using the 
algorithm briefiy described. 



5.5.2 Mapping to Langevin dynamics 

The goal of this subsection is to map our quantum field theory description of the rotors dynam- 
ics, which involves the two fields n*^^) and n*^^-* (see Sect. I3.3.4p . to an equivalent description in 
terms of Langevin dynamics. In the long-time limit of the coarsening dynamics, we establish 
that the equation of motion for the field n^^) is actually a Langevin equation driven by a colored 
noise $, the statistical characteristics of which are controlled by the self-energies of the fermion 
reservoirs. 

Let us take a step back and rewrite the effective action as it was before averaging over 
disorder. Making the assumption (we later check its consistency) that the Lagrange multipliers 
satisfy zf{t) = z^{t) = Zi{t) V i,t, the effective action reads 



n 



N 



N 



^E / dt f nr)(t) .nf (t) +iJ2j, n«(t) -nf (t) 

-\ I dt' Sf (t - t') nf )(t) • nf )(t') + i|dt' S|(t - t') n^\t) ■ n^\t') 
-i^,(t)n«(t)-nf)(t)} , (5.33) 



where introduced the real and symmetric matrix J' defined by J7ij = Jji/vN if j < i, J'ij = J'ji 
if j > i. Like the other components of this matrix, we set J'a to be taken from a Gaussian 
distribution with zero mean and variance J'^/N [we saw that the constraint nj(t) = 1 yields 
n^. (t) • n^. (t) = 0] . The total effective action adopts the quadratic form 

^S^K = -M^Jdt {inf(t)-^,(t) + i|dt'nf)(t).Sf(t-t')nf(t')} , (5.34) 



h 

where we introduced the N auxiliary fields ^f 

N 



7 = 1 -^ 



;df + z,it) ] 5i, - J^ 



6{t-t')-^^{t-t%Anf\t'). (5.35) 
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(2) 

By integrating over n^^ , we are left with 

N 

^5efr = -M^ IUtdt'Ut)-h:'f\t-t')Ut')- (5.36) 

From this Gaussian action, the quantity ^j(t) can be interpreted as a Gaussian random process 
with a zero average and variance (^i(t) -^jit'))^ = 6ijTi§{t — t') and eq. (15.351) as a set of coupled 
Langevin equations. This mapping is possible since the action of the rotor system, once the 
constraint on each rotor has been imposed through Zi(t) and Zi{t) is treated independently, is 
quadratic. In more general models the mapping is not exact, see e.g. the discussion in [52] . 

Under the further assumption Zi{t) = z{t), justified in the large M limit, the stochastic equa- 
tions (15.351) are rendered independent - apart from a residual coupling through the Lagrange 
multiplier - by a rotation onto the basis that diagonalizes the interaction matrix J'. Indeed, 
JT" being real and symmetric, it has A^ real eigenvalues Jo- with corresponding eigenvectors cr 
that constitute a complete and orthonormal basis of the space of rotor sites: cr • cr' = d^a' 
where • is the usual scalar product in this space. Let us collect all the rotors in the vector 
n = {n\ }je[i,Ar] and introduce its projections on the eigenvectors: n^ = n • cr. If we project 
eq. (15.351) onto cr, we are left with A^ uncoupled Langevin equations reading 



[^d^ - J. + zit)^ no(t) - j dt' Ei(t - t')no(t') = Ut) 



(5.37) 



with 

{Ut))i = , (Ut) ■ L'it'))^ = 5..' ^^{t - t') . (5.38) 

The noise statistics is peculiar because of the quantum origin of the environment: it has memory 
(colored) , and depends on T, eV, h. 

Two-time self correlation. Within the effective Langevin formalism, the two-time self cor- 
relation function defined in eq. (I3.23P reads 



C{t,t') = {n^it)-n^it')) , (5.39) 

where the average over disorder is realized by 

^T^-^= fdJ^pjiJ^) ■■■ , (5.40) 



and pj{Jcr) is the probability density of the eigenvalues of the interaction matrix J'. Following 
the analysis in [SO], the correlation function (I5.39P is expected to show a separation of time 
scales (at least in some parts of the phase diagram). This is usual in coarsening phenomena 
and corresponds to a stationary regime at short time- difference and an aging one at long time- 
difference with respect to a waiting-time dependent characteristic time. The stationary part of 

the correlation approaches a plateau at the Edwards- Anderson order parameter, q^A = (Ho-)? , 
that measures the fraction of frozen rotor fluctuations on time scales much smaller than this 
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characteristic time. The value of gsA depends on all parameters (T, eV, T, g). It is non- vanishing 
in the spontaneously symmetry-broken phase and continuously goes to on the critical surface. 
In certain cases it can be computed exactly. 

It is reasonable to expect that the long-time aging dynamics is determined by the low 
frequency (or long time) form of the Langevin equations only. The simplification arising in this 
asymptotic limit are discussed below. 

5.5.3 Long-time dynamics 

In the low-frequency, long time-difference limit, hw <^ T, the Keldysh self-energy can be ap- 
proximated by a constant [see, e.g., eq. fl4.20p in Sect. 14.2.31 for its exact expression in the 
ei? — )■ oo limit] 

Sf (r) ^ <5(r)Sf (a; = 0) > . (5.41) 

Similarly, we keep the leading contributions in the derivative expansion of E^: 

Sf (r) ^ Si(a; = 0)6{t) + r]6{T)dr , (5.42) 

with rj = (9t^;Im S|^(a; = 0) > 0. The Langevin equations read in this limit 

^d^n^it) + vdtn^it) = (j. - z{t) + S|(u; = 0)) n.(t) + |,(t) , (5.43) 

where t] plays the role of a friction coefficient and ^^rit) has white noise statistics: 

(Ut) ■ L'{t'))^ = S..'6{t - t') Sf (u; = 0) . (5.44) 

In the Langevin formalism, the kernel of an equilibrium white bath is given by the Einstein 
relation (known as the FDT of the second kind): {C,{t)C,{t'))^ = 2riT6{t — t'). Thus, the temper- 
ature T of the bath can be seen as the ratio of the diffusion coefficient of a particle embedded 
in that bath with the friction coefficient t] of the bath on the particle. For our reservoirs, 
in the low- frequency long time-difference limit, one can associate this ratio to an equivalent 
temperature T* 

T* ^ hm - , ^^^"1^, , , (5.45) 

the properties of which were discussed in Sect. 14.2.21 Thus, we confirm here that T* acts like a 
temperature in the sense that the effect of the (out of equilibrium) reservoirs on the long-time 
dynamics is the one of an equilibrium dissipative (Ohmic) bath at a temperature T*. This has 
been reported in different works and is at the root of the derivation of the stochastic Gilbert 
equation for a spin under bias [33j . 

We expect that as far as the long time dynamical behavior is concerned, the inertial term 
in eq. (I5.43P can also be dropped, thus leading to the equations: 

dtn^it) = X^{t) n,(t) + - |,(t) , (5.46) 

V 

where we introduced the shorthand notation Ao-(t) = [Jg- — Az{t)] /rj and ^z{t) = z(t)—T,^{u = 
0) and the spherical constraint is enforced by z{t). 
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This particular Langevin equation has been analyzed intensively in the study of the classical 
spherical Sherrington-Kirkpatrick model (or spherical p = 2 spin-glass model) and the results 
in [30] apply to our problem with T \-^ T*. The solution to eq. fl5.46p for a given disorder 
realization and noise history is 



n.(t) = n^(0) exp (j dr K{t)\ + - j dr ^,(r) exp (J dr' Kir'] 



(5.47) 



Copying results in [30j, the aging part of the correlation (in the limit t' ^ t — )■ oo) shows a 
simple aging scaling behavior 

C{t, t') ^ 2v^ gEA ^^l^n,-^,,, = Cit/t') . (5.48) 

The solution to eqs. (I5.46P leads to ^ea = 1 — T*{eV,T)/J. However, this result is obtained 
by taking the limit of relatively close times - with respect to t' - whereas, as we stressed, 
eq. (I5.46P is valid for the long time t' and long time-difference t — t' properties only. As a 
consequence, we expect the scaling result, eq. (15.480 . to hold at long times with the value of 
the Edwards- Anderson parameter not necessarily given by 1 — T*{eV,T)/J. Its computation 
requires a full solution of the equations of motion. 

We now focus on the aging dynamics in different parts of the phase diagram and argue 
that the Langevin dynamics of eq. (15.431) indeed provide a correct description of the dynamical 
evolution. 

Dynamics in the eV = plane. In this case, the Edwards- Anderson order parameter ^ea 
measures the static order parameter. The dynamic calculations based on the use of the quantum 
FDT to relate the correlation to the linear response in the stationary regime detailed in [29] , 
or the replica equilibrium computation in [35], can be easily extended to deal with a generic 
electronic bath in equilibrium. One confirms that ^ea = 1 at T = F = eV = and continuously 
approaches on the critical line Tc{T) for all values of g. The precise variation of gEA within 
the coarsening phase depends on the bath kernels. In the ei;' — )■ oo limit, the results in [22] 
apply also to our problem. The solution of the Schwinger-Dyson equations in the aging regime 
confirms that the scaling result, eq. (I5.48p . holds. 

Dynamics in the F = plane. Another interesting case is the effective overdamped 
Langevin limit obtained for F — )■ and (eV, T) in the coarsening phase. In this case drop- 
ping the inertial term in eq. (I5.43P is exact and not an approximation. 

Here the result gsA = 1 — T*{eV, T)/J can be shown to hold. The Edwards Anderson pa- 
rameter approaches one for T = V = T = and goes continuously to zero on the critical line, as 
in a second order phase transition. Consistently with the analysis of the critical surface derived 
from the QNESS phase (see Sect. I5.4.2p . one finds T*{Tc,eVc) = J. Numerical integration of 
the integro-differential equations of motion confirms that the scaling result, eq. (I5.48p . holds 
in the aging regime. 

Despite the fact that dropping the inertial term is exact, the equations (I5.46P are still not 
exact at all times. In particular, the initial conditions for this approximated equation of motion 
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Figure 10: (Color online.) Dynamics in the driven coarsening regime: numerical solution to Schwinger- 
Dyson eqs. (|5.6p and (j5.7p where the reservoirs have a half-filled semi-circle DOS with ep = IQJ- (a) The self 
correlation C{t,t') after a quench to 6* = 0.02.11 = 0.02, 7 = 0.2, 17 = 1 (in reduced quantities) shows first a 
stationary regime for short t — t' , then a slow aging regime where the time translational invariance is lost, (b) 
The self correlation C is plotted versus t/t' for two waiting times after two quenches into the coarsening region: 
9 — 0.02,1; = 0,7 = 0.2 and 9 = 0,v — 0.02,7 = 0.2. There is a double collapse of the curves. The collapse 
for the different t' proves the simple aging scaling C{t' /t) and the collapse for the two different quenches shows 
that T* ~ eV/2 plays the role of a temperature. The theoretical curve is the solution eq. (|5.48p with (/ea ~ 0.6. 

should be given by the state of the system a short while after the quench when the long-timescale 
description starts to be valid. Apparently, this delay seems to be not sufficient to significantly 
correlate the rotors with the interaction matrix J and, to any practical purpose no-(O) can still 
be considered "random" , at least as far as the Edwards- Anderson parameter is concerned. 

Dynamics in the T = plane. The zero-temperature plane is more difficult to deal with 
analytically. One is not entitled to use FDT since the system is driven by eV nor dropping the 
second time-derivative is exact. Furthermore, this is the case where the simplification leading to 
eq. (15.461) are more dangerous because of the power law tails appearing at T = in correlation 
and response functions. 

In order to check that the scaling result, eq. (I5.48p . holds we numerically integrate the full 
set of Schwinger-Dyson equations. 

In Fig. [To] (a) we show the decay of the two-time correlation function. For short time 
differences t — t' with respect to the waiting time t', there is a stationary regime depending 
on all control parameters in which the correlation approaches a plateau asymptotically in the 
time-difference. The plateau value is qea and measures the fraction of frozen rotor fluctuations 
on timescales much smaller than t' . Afterwards, there is an aging regime in which C depends 
on the two times explicitly. In Fig. [10] (b), we plot C against t/t' to prove that the simple aging 
scaling predicted analytically with eq. (I5.48P holds at these long times. Moreover, we show that 
the dynamics after a quench to^ = 0.2,f = are the same that the ones after a quench to 
9 = 0, V = 0.2, illustrating the fact that T* ~ eV/2 acts here like a temperature. 

Super-universality. It is remarkable that in the large M limit, the long-time dynamics of 
our model are exactly the ones of the classical fully connected p = 2 spherical spin glass. The 
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Figure 11: The integrated linear response, x(^j^') — Jt> dTR{t,T) against C{t,t'), for t' — 1024 and using t as 
a parameter. The curved part corresponds to the stationary and oscillatory regime with {t — t')/t' -^ while 
the straight line is for times in the monotonic aging decay of C. 

latter being a classical model in contact with an equilibrium bath (F = 0, eV = 0), the former 
being its quantum version in contact with a non-equilibrium bath (F ^ 0,eV ^ 0). The fact 
that the scaling functions are super-universal, in the sense that they do not depend on the 
external parameters T, eV, F once ^ea is extracted as a factor, can be understood as follows. 
First the fact that the non-equilibrium environment of our model gives rise to the same long- 
time dynamics than an equilibrium environment can be seen as a consequence of the Ohmic 
behavior of the reservoirs self-energy kernels at small frequencies (see Sect. I4.2.2p . Second, the 
fact that our quantum model shows a classical behavior at late times can be understood as a 
consequence of decoherence due to the dissipative (and Ohmic) bath. Furthermore, the effect of 
the temperature T on the long-time dynamics being irrelevant (in a RG sense) in the classical 
limit, one can expect the same to hold in the quantum case with respect to all parameters. 

This result has an interesting consequence. In the case of (large M) quantum 3d coarsening 
the classical-quantum mapping extends to space-time correlations and proves the existence 

-I /q 

of a growing coherence length R(t^) oc tw over which the rotors are oriented in the same 
direction. This real-space interpretation of aging unveils the connection with coarsening that 
was announced all along this manuscript. 

We found quite naturally that the long-time dynamics correspond to a Bose-Einstein-type 
condensation process of the A^ M-dimensional "vectors" n„ on the direction of the edge eigen- 
vector. The relaxation is controlled by the decay of p{Jcr) close to its edge. For Gaussian i.i.d. 
couplings p(Jcr) oc [(2^)^ — ■JaY^'^- This coincides with the distribution of the modulus of the 



Laplacian eigenvalues, p(fc^) ex {k 



2\d/2-l 



in c? = 3. For this reason all models with a square 



root singularity of the distribution of "masses" J^, as the ferromagnetic rotor model in (i = 3 
and the completely connected spin glass rotor model, are characterized by the same long-time 
dynamics. 
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5.5.4 Linear response 

It has already been noticed in Sect. I5.2l that the response function was somehow pecuhar since 
its equation of motion is decoupled from the one of the self correlation. Having argued that the 
long-time dynamics are governed by their classical counterparts, the linear response should also 
scale as in the classical limit. Therefore, the quantum fluctuation-dissipation relation between 
integrated linear response, xi^^t') — Jf^^" R(t,t") and self correlation C{t,t') approaches the 
classical one, x ~ ct -|- (gsA — C)/TeSi with an infinite effective temperature [36], T^s — )■ oo, 
as shown in Fig. [TTl The relations between integrated responses and correlation functions in 
other quantum problems that also approach classical-like form in the aging regime were shown 
in [SUEZ]. 

5.5.5 The Lagrange multiplier 

One should check the validity of a key assumption that was used to derive the phase diagram: 
the convergence of z{t) to an asymptotic value on the critical manifold. We first derive ana- 
lytically the asymptotic behavior (within our long-time approximation) of z{t) in the F = 
coarsening phase showing that this is indeed the case. Then we give numerical evidence that 
z{t) converges in the whole phase space. 

The condition C(t, t) = J dJ^ Pj{Ja) (ncr(i)-no-(t))5 = 1 reads after taking its time derivative 
and assuming furthermore that no-(O) is uncorrelated with a (no-(O) = no, V a), that is valid for 
random initial conditions (coming from infinite temperature for instance) 



= I dJ^pjM {dtn^{t)-n^{t))i: (5.49) 

= y"dJ.p,(J.)|n^A.(t)e2/od^Mr) + ^ I + 2Kit) J\t' e' 1'"^^" ^^(r") 1(5.50) 

Taking the derivative with respect to Hq yields 

= y"dJ.p^(J.)A.(t)e2/od^Mr)^ 



that can be recast into 



(5.51) 
Az{t) = I 5t In y" dJ^ pjM e''-'/^ . (5.52) 



Asymptotic behavior of z{t). By plugging the density of eigenvalues of an infinite (A^ — ?■ oo) 
and symmetric random matrix with Gaussian elements of variance J'^/N 




t-'rr 



Pj(J,) = — ^/l-(^) for J^ G [-2 J; +2 J] , (5.53) 



and zero elsewhere, we obtain 
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Figure 12: (Color online.) (a) Az{t) = z{t) - 
of the Jij matrix (here T = eV = T = Q.lJ,g 
and eV (dashed curve). 



(b) 



e, u 

'^%{'-^ — 0) quickly converges toward 2 J, the largest eigenvalue 
= 1, and ep — 10 J), (b) Dependence of z°° with T (plain curve) 



where Ji is the modified Bessel function of the first kind and first order, 
pre- asymptotic behavior for t 3> 77/ J 



We obtain, the 



z{t) ~2J + S 



UJ 



0) -r? — . 
' ' At 



(5.55) 



We just showed that inside the coarsening phase, the Lagrange multipher z{t) reaches an 
asymptotic value which is actually the critical value, z'^ = 2J + E§(a; = 0), calculated in 
Sect. 15.31 from the QNESS phase TTI equations without neglecting any term. The coherence 
between those two results somehow justifies the approximations made previously. In the ep — ;■ 
00 limit (reservoirs acting like an Ohmic bath) S/j(u; = 0) vanishes and we recover the same 
mechanism as in the classical case [30] . 

These analytical results are supported by the numerical analysis. Computed after the 
quench, the Lagrange multiplier z{t) quickly converges to an asymptotic value z°°. As an 
example, we plot in Fig. [T2] (a) the behavior of z{t) after a quench into the QNESS phase. The 
oscillations and the zero initial slope are signatures of the second and higher order derivatives 
in eq. (I5.37p . These terms were dropped in the analytical study of the long-time limit, see 
eq. (I5.46p . but the numerical integration does not neglect them. We give in Fig. [12] (b) the 
dependence of 2;°° with T and eV . It is quite clear that 2;°° is constant (and equal to z^) inside 
the critical surface and increases with T, F, and eV as soon as entering the QNESS phase. This 
justifies the assumptions made in Sect. 15.31 

To summarize the results, in the whole phase diagram z{t) always rapidly reaches an asymp- 
totic value z°° . Inside the QNESS phase, z°^ is a growing function of the parameters T, F, V^ 
whereas on the critical surface and inside the coarsening region, it is fixed to z^ . 
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Link between z{t) and the potential energy density One is interested in computing the 
energy density e{t) of the effective Brownian particle. It is given by 






(5.56) 



Using the solution fl5.47p for na-(t) at T* = 0, one has 

2e(t) = -n^ e-l /o^^^^^M f dJ„ J^p{J„) e^-^^*/'' . (5.57) 

By use of eq. (I5.52p . we obtain 

2e(t) = -|9i In f dJ^ pW e''^'/^ . (5.58) 

We recognize eq. (15.521) in the right-hand-side (rhs) of this last expression, giving finally 

e(t) = —Azit) . (5.59) 

This result is valid for any disorder density p{Ja)- For a non-zero T*, similar calculations give, 



see 



e(t) = l[T*- Az{t)] . (5.60) 



6 The current 



The physics of electric currents through mesoscopic quantum impurities in out-of-equilibrium 
settings has attracted a lot of attention in the recent years. The Kondo impurity is the canonical 
example of a strongly correlated system that has both been tackled experimentally [H] and 
theoretically by non-perturbative methods [32]. It is, to our knowledge, the first time that some 
fermionic reservoirs are coupled to a macroscopic disordered quantum system. In the previous 
Sections we analyzed the effects of the voltage drop on the system dynamics. In this Section we 
study the properties of the current that establishes between the two reservoirs. In particular 
we are interested in the possible influence of the rotors on the current. Is the current, that 
is rather easy to measure experimentally, able to give information about the dynamics of the 
rotors ? 

We recall the expression of the interaction Hamiltonian given in eq. (12.81) : 

^ N M Ns M 

HsB = -VM^J2T.m2< I'l'l^M <' i^myv +L^R]. (6.1) 

■^ i=l ij.=l k,k'=l l,l'=l 

From the point of view of the electric current, our model consists in two reservoirs coupled 
through time-dependent tunneling constants ^f(t). It is different from the usual quantum 
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impurity problems in the fact that the electrons cannot stay on the rotor system but only hop 
directly from one reservoir to the other. Furthermore, the quantum character of the system 
is not expected to play any significant role since its level spacings are smaller than any other 
energy scale in the large MN limit. The computation of the current will therefore lead to 
Landauer formula [13] a priori dependent on the rotors states. 

The electric current carried by the fermions flowing from the right to the left reservoirs is 

lR-.L{t) = -e (^) = ~{[H,Nr]) = -'^{[Hss.N,]) , (6.2) 

where — e is the electric charge of a fermion and Ni = Yliiki rnki'^Liki is the number operator 
of the left reservoirs. Hsb is the part of the total Hamiltonian H that couples the system and 
the reservoirs, see eq. (12. 8p . After straightforward algebra, we obtain 

iR-^dt) = -|(yM ^ Yl <'< [^'l^kl^R3k'v -L^r]). (6.3) 

ifikk'H 

In the Keldysh field theory formalism, this corresponds to the quantity 

In-.L{t) = \{l^_^^{t) + I-^_^,{t)) , (6.4) 

with 

iR^L^t) = -|(^^^ E ^"'<"(^) [r^kimiM^) -L^R]). (6.5) 

i^kk'll 

Expanding the action up to first order in the coupling constant g, we obtain an average over 
the rotors and the free fermions that are now uncoupled, that we note ( ■ • ■ )sb 

= ^^ (^) ' 5Z E E ^ / d^' <'<rn' {<\t)nf{t') (6.6) 

^ ab ifjbkk'lV juqq'mm' 
X [ruklitWR.k'vit) -L^R] WL,qmit')^'R,q>m'it') +L^R])sb. 

Averaging over the free fermions, we obtain 

iR^Lit) = -^MN{nw,f ^ 6 /" dt' ihG^\t, t') [ihGi;{t, t')ihG''^{t', t)-L^R] . (6.7) 

ab=± -^ 

G"-^ are the macroscopic Keldysh Green's functions for the rotors and G^^,^ are the Green's 
functions of the free fermions in the L/i?-reservoirs. This reads, after Keldysh rotations. 



lR-^L{t) = -^MN j dr C{t, t - r) nf (r) + R{t, t - r) Hf (r) 



(6. 
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with 






-2{nwcf Im 



^L^R ~^[^L^R +^L^r) 

^l'^r + '^l'^r \ ■ 



(6.9) 



The expression for the current given in eq. ( 16.8^ is quite generic. It is vahd as soon as the 
system and the fermionic leads are coupled with an interaction Hsb- The details of the system 
and the leads enter in the formula through their respective Green's functions. The formula was 
obtain after a first order expansion in the coupling constant g. The second order term like all 
the even order terms are zero by use of Wick's theorem. The third and higher odd order terms 
would have involved higher order correlation functions of the system. Plugging the expressions 
of the fermionic Green's functions G^ , G^, and G^ (a = L, R) that are given in Appendix lB.lt 
we get 



1 



u^iT) = ^inw^r{{ 



tanh(/3^^^i^)tanh(/3^^^^; 



n|(r) = -{hu^m 



tanh (/3 



2 



tanh (/3 



^R- I^R. 



sin|^^-^r))z.)^, (6.10) 



cos|^^^^r))i)«e(r), (6.11) 



h 



where the notation (( ■ ■ ■ )l)r stands for JJdede' pL{e)pji{e') ■■■ . One can check that the 
current vanishes when the bias voltage {eV = fiR — /i^) is set to zero. 



Linear conductance. We develop the current formula (16. 8p to the first order in eV and 
compute the hnear conductance 



IR^L 



(t)- 



-^MNeV [ dr C(t,t-r) 
^ Jo 



dn 



B 



[T 



\eV=0 



deV 



eV=0 



dnf 



deV 



, (6.12) 

eV=0 



One can derive for a flat half-filled DOS, p(e) oc Q^tp — |e — €f\), in the limit ei? — )■ oo (in that 
limit we expect the results to depend very little on the precise shape of the DOS) 



dnf (r) 



deV 
dnf (r) 



deV 



eV=0 



eV=0 



-7fg'^S{T) 



^ 2r 



Therefore the linear current very quickly goes from zero to 



'iJ-s-L 



it) 



2h 



MNg^ eV [n + h dr 



R{t,t-T) 



The dependence on the history of the two-time correlation function has disappeared 
second term in eq. (I6.15P goes to zero due to the rapid decay of the response function, 
the current quickly takes an asymptotic value 

Ir-^l = ^^MNg' eV . 



(6.13) 
(6.14) 

(6.15) 

and the 
Finally 

(6.16) 
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From this computation, it appears that the current only probes the very fast dynamics of the 
system it passes through and does not give information on the long-time dynamics. Since the 
short-time dynamics of the system are equilibrium ones even in the coarsening regime, the 
current cannot be used to tell in which regime the system is. An exact numerical integration 
of eq. fl6.8p supports these findings for other types of DOS, for finite values of ep and far from 
the linear regime. 

7 Conclusions and discussion 

In this paper we presented a detailed study of the quantum fully-connected rotor model driven 
out of equilibrium by a fermionic drive. We determined analytically the phase diagram of 
the model and we showed that a critical manifold, controlled by the value of the disorder 
strength, separates a QNESS with zero order parameter from an ordering phase with non-zero 
order parameter. We solved the equations that describe the dynamics in the different phases 
with a numerical integration and analytically by using various approximation schemes that give 
valuable physical insights. In particular, we showed that this (quasi) quadratic model maps to a 
set of Langevin equations with additive colored noise that describes the dynamics of the rotors. 
The nature of the noise is determined by the type of electron baths used and, in the driven 
case, the friction kernel and noise-noise correlation are not linked by any fluctuation-dissipation 
relation. By using this effective Langevin description we established the connection with the "id 
coarsening dynamics of the 0{M) model and we showed that the long-time ordering dynamics 
are in the class of the classical limit of our model without a drive, i.e., with the typical length 
growing as t^^'^ . 

Finally, we derived a generic expression for the current flowing through the system that 
involves a time-convolution between the characteristics of the system (through its correlation 
and linear response) and the ones the leads (through their retarded and Keldysh kernels). Inter- 
estingly enough, for the type of density of states used in the large ep limit the current depends 
only on the short-time difference (stationary) regime in which coarsening is not relevant. 

Future studies along these lines include the analysis of the fate of first order phase transi- 
tions, common in disordered quantum spin systems with multi-spin interactions when driven 
out of equilibrium. 
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Appendices 

A Conventions 

9 is the Heaviside step function. We choose 9(0) = 1/2, so that 9(a;) + 9(— x) = 1 V a; G M. 
We recall the identities 

/ '^^^xy ^ ^(^) ^^^ / ^^ ^^^) = 9(i/) , (A.l) 

where 5 is the Dirac delta function. In particular j_ dx 5{x) = 1/2. 

A.l Fourier transform 

The convention for the Fourier transform J-" that we use is 

/oo 
dr e+-- /(r) , 
-oo 

•^-H/M](r) = /(r) = /^^e~'^VM, 

The Fourier transform of the step function is 

jr[0(r)l(u;) = i pyi + 7r(5(u;) , (A.3) 

u 

where 'pv' denotes the principal value. Convolutions in real and Fourier spaces are defined by 



(A.2) 



{fog){r) ^ j dr'f{r')g{T-r') = T-'[{fg){uj)]{T), 
(/°^)M = f^fiuj')g{u-uj')=T[{fg)ir)]iu). 



(A.4) 



A.2 Heisenberg representation 

In the Heisenberg representation the operators evolve as 

Anit) = [/t(t)A(t)f/(t) . (A.5) 

with the unitary operator 

U{t)=Te-yo'it'Hit') ^ (A.6) 

and thus f/^(t) = Te"'^-'* "^^ ^i^), J and T are respectively the time and anti-time-ordering 
operators (see Appendix I A. 3 1) . For Hamiltonians H that do not explicitly depend on time we 
get 

Au{t) = e'^'/^A{t)e-'^'/^' . (A.7) 
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A. 3 Time-ordering operator 

On the real time axis, the time-ordering operator T rearranges operators with ascending times 
to the left: 

T Au{t)Bn{t') = An{t)Bu{t')Q{t - t') + C Bn{t')Au{t)Q{t' - t) , (A.8) 

with C = — 1 if both A and B are fermionic operators, ( = 1 otherwise. The anti-time-ordering 
operator T rearranges operators the other way round: 

t Anit)Bn{t') = An{t)Bnit')eit' - t) + ( Bn{t')Anit)e{t - t') , (A.9) 

On the Keldysh contour C, the position of an operator is specified by both the time and 
the branch index. By the notation Ai{{t,a), we denote the operator A in the Heisenberg 
representation at time t {t E [0,+oo[) on the branch a {a = ±). One can similarly define a 
time-ordering operator T^ that rearranges operators along the contour C represented in Fig. |5l 
The rules are 



(A.IO) 



JcAn{t,-)Bn{t',+) = An{t)Bu{t') , 

JcAn{t,+)Bn{t',-) = C Bn{t')Au{t) , 

TcAn{t,+)Bu{f,+) = An{t)Bu{t')e{t-t') + CBn{t')Au{t)e{t'-t), 

TcAnit,-)Bn{t',-) = Anit)Bu{t')Q{t'-t)+(Bnit')Au{t)Qit-t'). 

A. 4 Green's functions 

Let and 0^ be respectively annihilation and creation operators (bosonic or fermionic). In the 
field theory formalism of the Keldysh approach, we define the Green's functions as 

ihG''\t,t') = (0"(t)0^(t')) . (A.ll) 

a, 6 = ±, is either the complex conjugate (for bosons) or the Grassmannian conjugate (for 
fermions) of and the average is understood as 

(■■■) = y"p[0±,0±] ■■■ exp(^^S[<P^,^^]\ . (A.12) 

In the operator formalism the Green's function read 

mG''\t,t') = Tr [Tc Mt,a) (pl{t',b) ^h(0,±)] , (A.13) 

where 0h(^, cl) denotes the Heisenberg representation of the operator at time t on the a-branch 
of the Keldysh contour. Ph(0, ±) = ^(0) is the initial density matrix (normalized to be of unit 
trace) and its location on the + or — branch does not matter thanks to the cyclicity of the 
trace. Tc is the time-ordering operator acting with respect to the relative position of (t, a) and 
{t',b) on the Keldysh contour (see Appendix lA.Sp . 

One has, independently of the bosonicity or fermonicity of the field 

G''\t',t) = -G~^%t,t'y , (A.14) 

where the star indicates complex conjugate and a = —a. 
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B Fermionic bath 

We define the fermionic Keldysli Green's functions 

ihG'^\t,t')^{rm''it')), (B.i) 

where a, 6 = ±. Like for bosons [see eqs. fl3.2ip , one has 



G++{t,t') = G-+{t,t')Q{t-t') + G+-{t,t')Q{t'-t) 
G—lt,t') = G+-lt,t')e{t-t') + G-+lt,t')e{t' -t) 



(B.2) 



leading to the relation between Keldysh Green's functions 

G++ + G— = G+- + G-+ . (B.3) 

B.I Keldysh rotation 

We introduce the new fermionic fields 






(B.4) 



These definitions leads to 



ihG^^^\t, t') = {ip^^\t)f^\t')) = ih/A [G++ + G— + G-+ + G+-] = G 

ihG^^^^t, f) = {^^^\t)f^\t')) = i/2 [G++ - G— + G-+ - G+-] = -iG 

ihG^^^\t, f) = {ij^^\t)ipj^\t')) = i/2 [G++ - G— - G-+ + G+-] = iG ' 

i^G(22)(t,t') = (V'(2)(t)^(2)(t')) = i/fi [g'++ + G— - G-+ - G+-] = . 



K 
R 



(B.5) 



Where we defined, en passant, the Keldysh G^ , the retarded G'^ and the advanced G^ Green's 
functions in the same manner that we did for G and R in Sect. I3.3.4I Using relation (IB.3P we 
get 

G^ = ih/2 [G++ + G—] = ih/2 [G+- + G-+] , (B.6) 

G"^ = - [G++ - G+-] = [G+- - G-+] 0(r) , (B.7) 

G^ = [G++ - G-+] = [G+- - G-+] 0(-r) , (B.8) 



which are inverted as 



inG''^ = G^ + -(^aG'^-bG^) . (B.9) 



B.2 Symmetry properties under t ^ t' 

Using eq. (lA.14p . one establishes 

G'^ir) = -G^{-r)* , G^'ir) = G^'i-r)* . (B.IO) 

And hence in Fourier space 

G'^iio) = -G^ioo)* , G^H G M . (B.ll) 
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B.3 Free fermions 
B.3.1 Single free fermion 

The free fermion Hamiltonian is 

H = eil)^ . (B.12) 

Starting from the expression in terms of operators of the Keldysh Green's functions, 



',t (+1 



ihC^'it, f) = Tr [Tc i^uit, aWjj^t'. &)f?(0)J , (B.13) 

with a, 6 = ± and the grand-canonical density matrix ^(0) ex e~^^^^^^\ one computes 

ihG-+{e;r) = (l-n^)e-i^". ^ ' ^ 

np is the Fermi factor given by np{e) = (l + e^'-'^^^-*) . After the Keldysh rotation we get 



G"(e;r) = - tanh 1/3^1 e"^ 

G«(e;r) = le-^^-e(r) , (B.15) 

G\e;T) = le-t-e(-r). 

B.3. 2 Collection of free fermions 

For our left and right reservoirs, we consider continuous distribution (density of states) Pl{^) 
and pnie) of these free fermions. This yields to the Keldysh Green's functions 

G^\r) = ldep^{e)Gl\e-r), (B.16) 

with a = L,R. After a Keldysh rotation it yields 

G^{r) = /"de p(e)itanh[/3(e-/i)/2]e-^^^ = i(tanh[/3(e-/i)/2]e-^^0. , 

^^(r) = y"de p(e)-ie-i-0(r) = -^(e-i-),0(r), (B.17) 

G^{t) = Ide p(6)^e-i-e(-r) = ^(e-^-),0(-r), 

where we introduced a short-hand notation for the integration over energy levels. In terms of 
the Fourier transforms of p(e) it reads 

^^(r) = ^ 27rpir/h)eir) , G\t) = '- 27rp(r/n)0(-r) . (B.18) 
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B.3.3 Fourier transforms 



G^{u) = TT/itanh ( (3— — - ) p{huj) e R 



G^{u) + G^ico) = 2inp{huj) G iM . 
Since p(e) is real, one computes 

ImG^(u;) = 7Tp{hLo) . 



(B.19) 



(B.20) 



Thus we get, as a check, the grand-canonical fermionic fluctuation-dissipation theorem that is 
established generally in Sect. O 



G'^iuj) = ^tanh ( 13^^^^—^ \ Im G^'iuj) . 



(B.21) 



C Fluctuation-Dissipation Theorem 



In this Section we give a proof of the fluctuation-dissipation theorem both in its bosonic and 
fermionic versions. This theorem only holds in equilibrium and gives a relation between the 
Green's functions. In the grand-canonical ensemble, the initial density operator reads ^(0) oc 
exp{—P{H — fiN)), where A^ is the number operator commuting with H (in non-relativistic 
quantum mechanics), p is the chemical potential fixing the average number of particles. One 
can obtain the theorem for the canonical ensemble by formally setting p = 0. Let us consider 
a pair of either bosonic or fermionic operators, for instance creation and annihilation operators 
(j)^ and (p. Let us write the following Keldysh Green's function 



ihG+-{t,t')=Tv [Tc </'H(t,+)4(^',-)p(0) 
By resolving the time-ordering we get 

ihG+-{t,t') = ( TT\4>l{t')Mt)Q{0) 



id] 



(C.2) 



with C = +1 in the bosonic case and (^ = — 1 in the fermionic case. Using the analyticity of the 
Green's functions and then expanding 0h(^ + i/3^) = exp {—f3H) (piiit) exp (+/3if ), we get 



inG+-{t + iph,t') = CTr 

oc C Tr 



<(t') exp i-PH) Mt) exp {(3fiN) 



(C.3) 
(C.4) 



Since H and A^ commute and since for any operator f{N), one has (j)f{N) = f{N + 1)0, we 
have 



Mt) exp {/3pN) = exp {f3fi{N + 1)) Mt) 



(C.5) 
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and so 

inG+-{t + iph,t') = Cexp(/3/i) Tr Uji(t')^(0)(/)H(t) 

Using the cydicity of the trace, we come to 

ihG+~{t + iph,t') = Cexp(/3/i) Tr UH(t)0H(^')^(O) 
= C exp(/3/i) mG-^{t,t') . 



(C.6) 



(C.7) 
(C.8) 



If the system is in equihbrium, the tzme translational invariance of the previous equation gives 
the KMS relation: 

G+-{u) exp{l3huj) = Cexp(/3^) G-+{u) . 

Using eqs. (IB. 71) and fIB.Sp . we have on the one hand 



G^{u) + G^(w) = G+-(u;)(l - C exp{P{huj - /i)) . 
On the other hand eq. f lB.6P imphes 



(C.9) 
(CIO) 



ih 



G^{u) = -G+-(u;)[l + Cexp(/3(na;-/x))] 
These two last relations yield the grand-canonical quantum FDT: 

G'^iu) = h tanh (p^^^^^] Im G^fw) . 



(C.ll) 



(C.12) 



D Computing the self-energy 

D.l Derivation within the Schwinger-Keldysh formahsm 



In the Schwinger-Keldysh path- integral representation we had (see eq. (13.121) ) for the whole 
system (rotors and environment) 

Z[r7±] =y"2)[n±,V'±,V^^]e^^[-^'^^'^^l(n+(0),V'^(0)|f^(0)|n-(0),V'-(0)) , (D.l) 

At time t = 0, just after the quench, the initial density is assumed to be factorized: g{0) = 
Is ® gt^^O) O ^|°^(0) (see Sect. E2> yielding 

{n+iO),^l^''iOMO)\n-{0),^p-{0)) 

= 5(n+(0) - n-(0)) {i^tiOWnOM-Lm (V'i(0)|^r(0)IV'^(0)) • (D.2) 



The generating functional reads 



Z[rj 



±1 



V[n'^,n ]e« 



Ss[n+,n. ,r]] // gi^sBt^+.V" .V" :"■ -V" -V" ] \^\„ 
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(D.3) 



The index c' at the bottom of the integral is here to remind the constraints on the field inte- 

9 9 

gration, namely n.f{t) = n^it) = 1 and n^(0) = n~(0) V i. We introduced the average over 
the free environment composed of the two reservoirs: 

{{■■■)l)r^ jviijj^A^] ■■■ ei^^ei^H 

We now develop the coupling e^'^^^ up to the second order, 

(( e^^- ),)n ^ 1 + ^(( SsB )l)r - ^(( Sl^ ),)n . (D.5) 

The first order term is zero. The second order term reads 

^ '^ ^ ab=± "^ •^'^ ij=l kk'qq'=l lJ.u=l U'mm'=l 

X (( [ruMmmk'i'it) + Li^R] [^Pi,,^{t')^%,Ut') + L ^ R] )l)r . (D.6) 
Developing the term on the second line, we obtain 

(( [^hkiit)rmk'i'it) + L ^ R] m,,mit')^%,Ut') + L ^ R] )^)« 

= (( rR^Mit)rukn'minn^(t')^Rn'rn'i^ +L^R ) l) R 
= -{i'hk'At)i'lgm{t'))L {i^Rjq'm'it')rRikMR + L^R 

= S.jh'MSi'mSim'h' [GfAt, t')G%{t', t) + L^R] . (D.7) 

With the free fermionic Green's functions defined on the Keldysh contour as i^G^^(t,t') = 
{4'ki't)'ipki't')) a for a = L,R, a,b = ± and where k labels the electron's energy. Expression 
(ID. 61) now reads 



/ f. \2 n i^oo N Ns M M 

( Si, ),n = ^'mI^) Y^ah dtdt' Y.Y.Y1 E<"(^)<'(^') -«'^^'^ 

^ *^ ah=± •' -^^ i=l kk'=lfiu=lU'=l 

x[Gf,,{t,t')G''^,{t',t) + L^R] . (D.8) 

By using the property Tr cr^cr'^ = 5fj,u, we get 

2 . .oo N 



{ Si, ),n = Mh' (^) Yl «^ //°°dtdt' X:<(t) ■ nlit') 

X Y. [^^'fc'(^' t')G'Skit', t) + L^R] . (D.9) 

kk' 

Finally expression f lD.5|) can be recast into 

(( e^^«^ )l)r ^ e^^'sl , (D.IO) 
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with 



-I p r+oo 

Sfl[n\n-] ^ —M J2 // dtdt' S^^(t, t') J]<(t) ■ n,^(t') , (D.ll) 

where the exponent (2) is here to recall that we developed until second order and with the 
self-energy 

E^^(t,t') = -abih {hw^Y [Gfit,t')G'^%t',t) + G'i{t,t')G'i^it',t)] , (D.12) 

where the Keldysh Green's functions of the fermions in the a-reservoir {a = L, R) are given by 

Gl\t, t') = y"de„ p^leJG^'le,; t - t') = G^^(t - t') . (D.13) 

Pa(e) is the density of states in a-reservoir and G'^{e; r) are the Keldysh Green's functions of 
a free fermion with energy e in equilibrium in the a-reservoir (see Appendix IB. 3. ip : 

ihG+-{e;T) = -n„(e)e-i^- , 

ihG-+{e;T) = [l-n„(e)]e-i^^, /j^ ^^x 

inG'++(e;r) = i;iG'-+(e; r)e(r) + i;iG'+-(e; r)0(-r) , ^ ' ^ 

ihG--{e; t) = ihG+-{e; t)&{t) + mG-+{e; T)e{-T) , 

with the Fermi factor na{^) = (1 + e'^"''^'^"'^'"-')"^. It is clear then that the self-energy is time 
translational invariant: Il'^{t,t') = T,'^{t) with t = t — t'. Moreover S^''(r) is a symmetric 
matrix with respect to time and Keldysh indices: 

Eg'(r) = S^"(-r) , (D.15) 

Using the time reversal property eq. (lA.14p of the Keldysh Green's functions one also establishes 

S^^(r)* = -Sf (r) , (D.16) 

where we note a = —a. 

After a Keldysh rotation of the rotors coordinates, it yields 

• -j (. poo ^ 

-4^][n«,n(^)] = -M Y. li dtdt'S:e^(t,t')$^nKtK(t'), (D.17) 

with 



^ ^ J Jn 



(D.18) 



s^,^^) = -i/h [s^+ + s+- + s^+ + s^-] = . 

which is inverted as 

ihE'S = -ahUf^ - ^ (aHf^ + 6S(;^)) . (D.19) 
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D.2 FDT check 

We checked that the ferinion-reservoir self-energy satisfies the bosonic FDT. This is only valid 
when the reservoirs constitute an equilibrium bath, i.e., (3^ = /3r = /3 and fiL = f^R = /^o 
{V = 0). Note that distribution functions piiyj) and Pr{uj) can be different although the proof 
given below uses piifjj) = Pr{^) = p(e) for simplicity reasons. The goal is to check 

Sf (c) = h coth (f3^) Im Si(a;) = /Icoth (/3^) ^^^ ^ ^^^ ^^^ . (D.20) 

We first develop the term in the LHS, then we do the same with the RHS to prove their equality. 

SfH = TFS|(r) 

= -2{nuj,f TF {G^G^* - h^A [G^G^* + G^G^*] } 

= -2{niu,y TF {G^G^* - h^A [G^ + G^] [G^* + G^*] } , (D.21) 

where we used the nullity of cross terms of the type G^G^ since G^ oc 0(t) and G^ oc 0(— r). 

Sf (w) = -2{hw,f {G^ o G^'* - fi^/A [G^ + G^] o [G^* + G^*] ] , (D.22) 

where o is the symbol for the convolution (see Appendix |X]) and G^*{u) stands for the Fourier 
transform of G^{t)*. Since we easily obtain 



G^{u) + G^{uj) = 2i7rp{hu) , 
G^*(a;) + G^*(u;) = -2inp{-huj) 



(D.23) 



and 



G^M = vrMM tanh (/3^1™) , 
G^*{u) = nhp{-hu) tanh {P^^^f^) , ^ ' 

we get by replacing in ( JD.22J) 

j:^{io) = -2{nLo,f{7rh)^ 

X {[p(na;)tanh(/3^^)]o[p(-;k^)tanh(/3^^^)] - [p{hu)]o[p{-huj)]} 
= -2{nL0cf{ixh)j%p{t')p{t' - hio) jtanh (/3^) tanh f/3'—Il^p^\ - l\ (D.25) 

= -Tih{huJcf coth {I3^f) /de' p(e')p(e' - hw) jtanh (b ''-^^-^A - tanh ((3'- 

where we used the trigonometry relation 

tanh X — tanh y 

tanh [x-y) = — . 

1 — tanh X tanh y 

Let's now calculate the RHS of (ID.20p . 

[Si + S^] (c) 



2i 



i(nwc)^ TF {G^G^* + G^G^* + G^G^* + G^G^*} 

i{huj,f TF {(G^ + G^)G^* + G^(G^* + G^*)} 

i{nw,f { [G^ + G^] o [G^*] + [G^] o [G^* + G^*] } , (D.26) 
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giving 

X {[p(M] ° [p(-Mtanh(/3^^)] - [p(^)tanh(/3^)] o [pi-hw)]} 
= -7rh{hu,f{27rh) coth (/3^) / de' p(e')p(e' - hu) jtanh (/j ^"^^"^" ^ - tanh (/3^) } . 

We recognize here the development (1D.25I) of S^. We just proved that the bosonic FDT is 
satisfied provided that the two fermionic reservoirs have the same temperature and chemical 
potential. They can have a different density of states. 

E Dynamics 

E.l Quadratic effective action 

One can render the effective action quadratic at the price of introducing new fields. For a given 
i and a given pair of (r, /i, t) and (s, u, t'), the identity 

1 = JdQ/,t{t,t') S {nr{t)nr{t') - Q:,t{t.t')) , (E.l) 

becomes, after using the integral representation of the delta distribution (see Appendix lAl). 



loc I dQ,;t{t,t')dX,;t{t,t') exp(-i-X,;t{t,t'){nr{t)nr{t')-Q,;t{t,t'))) . (E.2) 



,M 

y 

Introducing similar identities for all possible pairs of {r,fi,t) and (s,z/, t'), we obtain a path 
integral over twqj fields <5j^t(t,t') and Aj^^(t, t') that are symmetric in the Keldysh indices, 
times and rotor components: Q^f,^^{t' ,t) = QiJZ(t,t') and \t,^^{t' ,t) = \l,l{t,t'). The effective 
action is now also a functional of Q and A and reads 



\s^^ = -y E jh^^^' EE^'w [o<:(t,t') +iA/:(t,t')] <^(t') (e.s) 

r,s=(l),(2) "^ "^ i fif 



2N 



a i rs i ^v 

+ boundary terms , 



^There are N {AP K'^ + M K) /2 of each of these fields, where K = 2[s the number of possible Keldysh indices. 
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where we introduced the operator Opj^ll(t,t') defined as 



^dj + W'^(t) 



2... 



W^,Sf(t',t) 



a=± 



op.l!'\t,t') ^ ^5,jit-t')j2<it). 



a=± 



Opi^l{t,t') is symmetric in the Keldysh indices, times and rotor components: Op^l^'^it' ,t) = 
OPial{ti t')- The functional integration over n^^ is now quadratic and can be performed, leading 
to 

^5eff = -^Tr lnM(Op + iA) (E.5) 

n 2 



2 T\/r2 



rM 






2N 

a i rs i ^u 

where the trace in the first term is spanning the whole space of indices, namely rotor sites, 
Keldysh indices, times and rotor components. 

E.2 Saddle-point evaluation 

In this subsection, we evaluate in the limit M — )■ oo the saddle-point equations with respect 
to the dummy fields we introduced previously, namely Aj^*(t,t'), QilZi t,t') and z1(t). The 
fluctuations around the saddle are neglected. In particular, using eq. flE.ip we have the identity 
(see the definition of Green's functions in Sect. I3.3.4p 

Q,;t{t,t') = ihG,l^^{t,t'). (E.6) 

Along the lines we prove that the solution in the saddle is 0{NM), like the starting Hamiltonian. 
The saddle-point with respect to Xi^lit.t') yields 

^^''^ -\ Tr ,, J,^^, lnM(Op + iA) + i^Q-,(t,tO = , (E.7) 



<5A,-(t,t') 2 5A,-(t,t') 

giving in matrix notations 

\Op + i\)-^ = MQ , (E.8) 

where the symbol * represents the transposition. Since all operators in the last equation are 
symmetric by definition, we get 

Op + iA = ^g-^ (E.9) 
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The saddle-point equation with respect to Qilil(t,t') yields 



^Ktit,t') 



N 



Y.Q^^^^) ^ 



(E.IO) 



where (2) = (1) and (1) = (2). The RHS of this last equation being site-independent, Aj^^(t, t') 
does not depend on i: Xijll{t,t') = X'^^it.t'). Equations (IE.9P and (lE.lOp imply 



opr + 



J'^M 



Eq, 



j 



1 

M 



Q 



-l rs 



. 



(E.ll) 



The saddle-point equation with respect to z!^(t) yields to the two equations: 



fi- 



ll ipp+^>r'):: (M) + 'i ([o.+iA]-)^r; (m) = m . 



U (22) 
ififi 



(E.12) 



This is nothing more than the constraint that rotors should have a unit length. However, A 
being site-independent, it is clear from these equations that it has to be the same for Op. Finally 
at the saddle. Op, Q and z are site-independent (homogeneous) so we can get rid of the sites 
indices: Op,lt{t,t') = Op^;,{t,t'), Q,;^(t,t') = Ql,l{t,t') and zfit) = z%t). Equation dEIID 
becomes 

Op''' + .PMQ^' - ^Q'^'' = . (E.13) 



Since from its definition (]E.4p Op^^^{t,t') oc 5fj_y, the previous equation tells us that it has to 
be the same for Q[^'^^{t,t') so we can get rid of all the rotor component indices. Multiplying by 
Q'^^(t',t"), and summing over s and t', we get 

f dt' J2 Op'^t. t')Q'%t', t") + J^MQ'\t, t')Q'%t', t") - ^6rJ{t - t") = . (E.14) 

The macroscopic Green's function reading ihG'''{t,t') = MQ'''{t,t') we obtain 

f dt' Yl Op''{t, t')ihG'''{t', t") + JhhG'\t, t')ihG'''{t', t") - 5rJ{t - t") = . (E.15) 

E.3 Schwinger-Dyson equations 



The (r = (2),v = (1)) component of eq. ( lE.lSp gives a complex equation the real part of which 
yields 



z+{t) = z-{t) = z{t) Vt, 
and the imaginary part of which is the dynamic equation for the self-correlation: 
1 d^ 



Vdt 



2+^W 



C{t,t') 



dt"j:^{t,t")R{t',t")+ / dt"E^{t,t")Cit",t') 

Jo 



(E.16) 



(E.17) 
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where we introduced 

S^ = j2C + Ef, S^ = j2i? + E|. (E.18) 

Similarly, the (r = {2),v = (2)) component of eq. f lE.lSP yields the equation of motion for the 
self-response: 



1 92 



R{t, t') = 5{t - t') + / dt" S^(t, t")R{t'\ t') . (E.19) 

Jt' 



The (r = (l),f = (1)) component of eq. ( lE.lSp leads to the same equation and the (r = (1), f = 
(2)) component expresses = 0. Setting t' = t in eq. ( lE.lTp we obtain the expression for the 
Lagrange multiplier 



Equations (IE.17P and (IE.19P together with eq. flE.20P constitute the Schwinger-Dyson equations 



that fully determine the dynamics of the interacting system. 
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